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Infrared  Fluorescence  of  Electron 
Irradiated  CO2  in  the 

Presence  of  N2.  Ar,  and  He  j 


I.  IMIIOIMI nON 

The  I.ABCEDK  facilit-  at  the  Air  Force  Geophysics  L.ahoratory  (AFGF)  is  used 
to  observe  the  fluorescence  beiiavior  of  electron  irradiatert  gases.  In  this  experi- 
ment a constant  pressure,  continuously  flowing  gas  is  irradiated  by  a pulsed,  col- 
limated electron  beam  of  known  current  and  V'oltage.  As  the  electron  beam  passes 
through  the  gas,  th.e  beam  electrons  undergo  collisions  with  the  gas  molecules 
which  c sult  in  beam  energy  loss  and  spreading.  Some  fraction  of  the  energj  lost 
by  the  ( lectron  beam  will,  eitiier  directly  or  indirectly,  appear  as  fluorescence  in 
the  gas.  In  the  l.AHi  KDK  experimeiV,  thi  . fluorcsoence  is  observed  through  a 
window  with  the  optical  axis  of  the  detection  system  at  right  angles  to  the  beam 
axis. 

I'Or  the  ca-.e  of  infrared  fluorescence,  the  primary  diagnostic  is  a Micholson 
interferometer  whose  signal  is  channeled  tjuxiugh  a modified  lock-in  amplifier. 

The  interferometer  in  operated  in  a mode  prcividing  time  resolved  spectra,  as 
previously  described,  “ thus,  hign  spectral  resolution  infraied  histories  of  (he 
fluore  icing  gas,  both  during  and  after  electron  beam  excitidion,  can  be  developed. 

• (Heci  ived  for  publication  (!  September  li'VVl 

1.  Cook,  F.ll.,  and  Murpliv,  II.  K.  <10761  A Synchronous  Sigtial  Processing 
1 cchnidue  for  llepelitive  Arbitrary  Waveforms,  AFC'lIl  ,-TH -76 -00:i  5. 

• 2 Murphy,  II.  F.,  C'.iok,  !■'.  II.,  and  Sakai,  11.  (1'175)J.  Gpt.  Soc.  Amcr,  6a:li00. 
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Prelimitiarv  <ibdefvations  of  this  type  in  electron-i fradiated  gas  mixtures  of 

extiibited  spectra  typical  of  highly  vibrationally  excited  CO,  a species  not 
present  in  the  ambient  gas  mixture.  The  remainder  of  this  text  is  devoted  to  an 
analysis  of  tliese  results. 


2.  IMII  I IMINMO  lt\T\  VWI.^SIS 

Preliminary  measurements  of  the  infrared  fluorescence  resulting  from  the 
electron  irradiation  of  N2/CO2  mixtures  were  performed  utilizing  a two  liter,  gold- 
coated,  integrating  sphere  in  order  to  maximize  the  output  intensity.  Typical  con- 
ditions for  the  measurements  were  Pn2  torr,  PCO2  torr, 

beam  voltage  V .32-50  kV,  beam  current  I - 0.  ()-3  mA,  beam  pulse  time 
Tp  = 2.  5 - 5.  0 msec,  and  a beam  duty  cycle  of  ^ 20  percent.  A typical  measure- 
ment of  the  observed  time  dependent  fluorescence  is  shown  in  Figures  1 and  2.  In 
Figure  1,  the  observed  time  dependent  infrared  fluorescence  of  a typical  case  is 
shown  over  the  wavenumber  range  of  1700-2500  cm  \ The  spectral  resolution  of 
this  data  is  =10  cm  Note  that  the  data  shown  is  raw  detector  signal  and  has  not 
been  corrected  for  detector  response  and  system  transmission  losses.  The  two 
dominant  spectral  features  in  this  wavenumber  region  are  the  002(12)  band  between 
2300-2400  cm  \ and  what  is  apparently  fluorescence  from  the  fundamental  band 
of  highly  vibrationally  excited  CO  in  the  wavenumber  range  of  1750-2200  cm  In 
the  particular  case  shown,  the  electron  beam  pulse  duration  was  4.  fi  msec  and  it 
can  be  seen  that  the  002(1^2)  intensity  continues  to  rise  in  that  time  span  whereas 
that  of  Co  appears  to  approach  "steady  state"  within  1-2  msec.  After  beam  termi- 
nation the  fluorescence  decays,  the  CO  decays  more  rapidly  than  the  CO.,.  This 
can  be  seen  best  in  Figure  2 where  the  same  data  is  plotted  in  reverse  time  in 
order  to  more  prominently  exhibit  the  decay  behavior. 

In  general,  the  data  shown  in  Figures  1 and  2 are  typical  of  observations  made 
over  the  range  of  conditions  considered.  In  particular  for  all  cases  (a)  the  CX) 
intensity  is  of  the  same  order  as  that  of  the  C02(>'3)  b^ind  even  though  CO2  pres- 
ent in  the  ambient  gas  and  CO  i.s  not,  (b)  the  CO  fluorescence  reaches  steady  state 
more  rapidly  than  CO2,  and  (c)  the  peak  CO  intensity  is  independent  of  CO,,  concen- 
tration. Furthermore,  in  all  cases  the  CO  is  highly  vibrationally  excited.  It  can 
be  shown  that  for  a Boltzmann  distribution,  the  vibrational  temperature  T mav  be 

V * 

deduced  from  the  observed  spectra  by  the  relationship 
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Figure  1.  Time  Dependent  Fluorescence  in  an  N2/CO2  Mixture.  Pressure  listed 
in  torr.  Beam  conditions  v = 32  kV,  I = 0.  6 mA 


Figure  2.  Fluorescence  in  an  N2/CO2  Mixture  vs  Hover 
msec.  Conditions  as  in  Figure  1 
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where  T id  the  characteridtic  vibrational  temperature,  = 3100‘’k  for  CC),  and 

V id  the  vibrational  level  at  which  the  intendity  id  highedt.  For  the  dpectra 
max 

dhown  in  Figured  1 and  2.  V = and  thus,  if  the  CO  vibrational  didtribution 
^ max 

were  lloltzmann,  the  obderved  spectra  would  corredpond  to  a vibrational  tempera- 
ture of  = 19,  000“k.  Thid  id  to  be  compared  with  the  gas  trandlational  temperature 
which  id  nominally  ambient. 

In  principal  data  such  as  that  shown  in  Figures  1 and  2,  it  can  provide  a sig- 
nificant amount  of  kinetic  information.  In  particular,  rate  constants  for  the  vibra- 


tional  exchange  processes. 

CO(v)  + N2(v  = 0)  - CO(v  - 1)  4 N2  (v  1) 

(2) 

CO(v)  + CO2  -CO(v  - 1)  + 002(1/3) 

(3) 

could  be  determined  for  a wide  range  of  CO  vibrational  levels  from  the  CO  fluores- 
cence decay  observed  after  beam  termination.  Such  analyses  would  unfortunately 
be  complicated  if  the  CO  vibrational  densities  were  sufficiently  high  so  that  intra- 
mode  vibrational  exchange  processes,  that  is, 

CO(v)  + CO(v')  . CO(v  - 1)  + CO(v'  + 1)  , (4) 

played  a role  in  the  relaxation  processes.  In  any  event  once  the  manifold  of  decay 
rates  were  determined,  the  CO  fluorescence  observed  during  electron  irradiation 
could  be  analyzed  in  order  to  specify  the  vibrational  level  dependent  production/ 
excitation  rates  for  CO(v).  Furthermore,  for  the  cases  considered,  the  C02(i'3) 
mode  will  be  equilibrated  with  N2  (v  = 1)  through  the  rapid  intermode  vibrational 
exchange  reaction 

N2(v  = 1)  + CO2  <^002(1/3)  + N2(v  = 0)  (5) 

which  is  nearly  thermo-neutral.  Thus  if  the  CO^iv^)  vibrational  temperature  can 
be  determined,  the  nitrogen  vibrational  temperature  will  also  be  specified. 

Nitrogen  is  vibrationally  excited  directly  by  beam  electrons  (and  secondaries) 

e + N2(v  =0)  -►  e + N2(v  = 1)  (6) 

and  thus  a determination  of  the  N2  vibrational  temperature  could  lead  to  a determi- 
nation of  the  global  rate  constant  for  reaction  (fi).  In  particular,  it  would  be  desir- 
able to  determine  what  fraction  of  the  beam  energy  is  converted  to  vibrationally 
excited  nitrogen.  (As  will  be  shown  later,  such  an  analysis  is  complicated  by  the 
fact  that  the  C(l2  radiation  within  the  test  chamber  is  optically  thick). 
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'1  he  kinetic  interpretation  of  this  data  will  of  course  be  complicated  by  the 
fact  that  there  are  three  vibrationally  excited  species  present,  CO,  and  CO^, 

which  may  interact  with  each  other.  Therefore,  in  order  to  simplify  the  vibra- 
tional phenomenology  as  well  as  to  more  readily  elucidate  the  CO  production 
ineclianism,  fluorescence  measurements  were  performed  in  electron  irradiated 
mixtures  of  Ar/CO,,  and  He/CO,,.  These  measurements  were  performed  for  the 
same  beam  conditions,  pressures  and  CO,,  diluency  as  the  measurement  in  the 
N,,  CO,,  mixture.  A comparison  of  tlie  observed  time  dependent  spectra  in  the 
three  cases  is  shown  in  Figure  3 while  a comparison  of  the  spectra  at  one  time  is 


Intensity  Decay  of  CO, 
with  Various  Gases 


Figure  3.  Fluorescence  in  X/CO-j  Mixtures  vs  Heverse  Time  in  Steps 
of  0.  1 msec,  Pv  30  torr,  P,.,,*”  0.  ,3  torr.  Beam  conditions  as  in 

T^-  t to,. 

Figure  1 3 
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1 shown  in  Figure  4.  As  can  be  seen,  highly  vibrationally  excited  CO  is  observed 

in  all  three  cases,  although  the  CO  vibrational  population  distributions  vary  from 
case  to  case.  Note  that  the  C02(i'3)  radiation  decreases  more  rapidly  in  the  He 
and  Ar  mixtures  than  in  N2.  This  is  because  the  vibrationally  excited  nitrogen 
acts  as  an  excitation  source  for  002(1’ .j),  thus  effectively  increasing  its  apparent 
lifetime. 


I The  observations  in  Ar  and  He  are  quite  useful.  In  particular,  they  demon- 

' strate  that  the  production  of  vibrationally  excited  CO  is  not  related  to  the  presence 

of  N2.  Furthermore,  the  rate  constants  for  vibrational  relaxation  of  CO  by  CO2 
may  be  more  readily  deduced  from  the  Ar/C02  and  He/C02  data.  These  results 
could  then  be  used  in  interpreting  the  N2/CO2  observations. 

Before  any  such  analyses  can  be  performed,  however,  the  kinetic  behavior 
occurring  within  the  test  chamber  while  the  beam  is  on  must  be  diagnosed.  As  the 
electron  beam  propagates  through  the  gas  it  loses  energy,  slows  down  and 
spreads.  The  degree  of  spreading  and  energy  loss  is  determined  by  the  thickness, 
that  is,  gas  number  density  times  length,  of  the  gas  target.  The  lost  beam  energy 
is  transferred  to  the  gas  causing  molecular  excitation,  dissociation,  ionization, 

I and  heating.  Through  such  phenomena,  species  not  originally  present  in  the  gas 

are  created  and  these  can  effect,  as  well  as  produce,  the  resulting  fluorescence 
histories.  Furthermore,  since  the  beam  spreads  and  slows  down,  the  electron 
j excitation  rates  are  spatially  nonuniform  whereas  the  fluorescence  observations 

i are  integrated  over  a two-liter  volume.  Lastly,  species  created  by  the  electron- 

I beam  can  interact  with  each  other,  or  with  ambient  species,  to  produce  additional 

products.  The  electron/ion  kinetics  of  a typical  case  is  considered  below  in  order 
to  determine  typical  electron  beam  species  production  rates  and  the  ultimate, 
gaseous  products  resulting  from  electron  irradiation  of  N^/C02  mixtures. 

I 2.1  (')lertrnn 'Molecule  Keac'tions 

Neglecting  excitation  processes,  the  dominant  electron-neutral  reactions  in 
irradiated  N2/CO2  mixtures  are  as  follows: 


(1)  Ionization  of  N2 

e + N2  2e  + N2^ 

(7) 

- 2e  + N + N"^ 

(8) 
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(2)  Ionization  of  CO2 

e + CO2  -*  2e  + C02^ 

(9) 

2e  + CO  + 0^ 

(10) 

- 2e  + 0 + CO^ 

(11) 

2e  + 20  + C”^ 

(12) 

(3)  Dissociation  of  N2 

e + N2  -»e  + N-i-N 

(13) 

(4)  Dissociation  of  CO2  j 

e + CO2  - e + CO  + O <14) 

(5)  Etc  . . . 

The  cross  sections  for  the  first  three  of  these  processes  are  shown  as  a function 
of  electron  energy  in  Figures  5-7.  The  data  of  Figure  5 are  taken  from  the  work 
of  Rapp  and  co-workers,  ^ and  that  of  Figure  6 from  Crowe  and  McConkey.  ^ 

The  data  of  Figure  7 represents  the  difference  between  the  total  N2  dissociation 
cross  section  reported  by  Winters*’  and  the  dissociative  ionization  cross  section  of 
Rapp  et  al.”^  Little  information  is  available  on  the  direct  dissociation  of  C02. 

This  process  has  been  dropped  from  the  analysis  inasmuch  as  the  mixture  ratio 
under  consideration  is  N2/CO2  of  100/1,  and  the  CO2  dissociation  cross  section  is 
anticipated  to  be  less  than  that  for  ionization  except  at  the  lowest  electron  energies. 

It  is  difficult  to  use  these  cross  sections  directly  to  calculate  the  various 
species  production  rates  inasmuch  as  such  production  occurs  through  an  electron 
cascade,  and  thus  the  effective  rate  involves  an  integral  over  the  electron  energy 
distribution.  The  average  electron  energy  required  to  create  an  ion  pair,  has, 
however,  been  measured  for  a number  of  gases  including  N2;  thus,  the  fractional 

3.  Rapp,  D. , and  Englander-Golden,  P.  (191)5)  J.  Chem.  Phys.  43:14(14. 

4.  Rapp,  D. , Englander -Golden,  P.  and  Briglia,  D.  (1905)  J.  Chem.  Phys. 

«:4081. 

5.  Crowe,  A,,  and  McConkey,  J.W.  ( 1974)  J.  Phys.  B7:349. 

0.  Winters,  H.  F.  (1900)  J.  Chem.  Phys.  44:1472. 
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Figure  5.  Ionization  Cross  Sections  for  N2 


Figure  (i.  Ionization  Cross  Sections  for  CO, 
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Figure  7.  Dissociation  Cross  Section  for  N2 

species  created  per  ion  pair  can  be  estimated  by  ratioing  the  cross  sections  of  the 
individual  process  to  that  for  N2  ionization.  This  has  been  done  at  an  electron 
energy  of  100  eV,  which  is  well  above  the  threshold  energy  for  the  reactions  con- 
sidered. The  results  of  this  analysis  are  presented  in  Table  1 for  the  case  of  an 
N2/CO2  ratio  of  100.  As  would  be  anticipated,  the  dominant  ions  are  N2^  and  N^; 
however,  several  other  trace  ions  are  created  and  these  could  play  an  important 
role  in  the  further  analysis. 

2.2  Ion  MuUm'uIi'  llrarliiins 

The  ions  created  by  direct  electron/molecule  collisions  will  ultimately  recom- 
bine with  slow  secondary  electrons  via  rapid  two  body  mechanisms  such  as 

AB*  + e * A + B . (If)) 

“7 

Such  recombination  reactions  typically  have  rate  constants  of  order  10  cm'  /sec 
at  thermal  energies.  However,  since  the  electron/ion  concentrations  in  the  gas 
will  generally  be  less  than  10  /cm  , the  characteristic  time  for  ion  destruction 
will  be  greater  than  10  psec.  As  will  be  shown,  a number  of  ion/molecule  reactions 
will  occur  in  the  gas  on  time  scales  much  shorter  than  10  ;jsec.  The  net  effect  of 
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Table  1.  Ion  Production  Per  Ion  Pair  N2/CO2  = 100 


Spec  ies 

Fraction/Ion  Pair 

Additional  Reaction 
Products 

^2" 

-0.75 

- 

n" 

-0.24 

N 

C02^ 

-0.01 

- 

0" 

0.0017 

CO 

0 

0 

+ 

0.0010 

0 

0,  00075 

20 

these  reactions  is  that  species  with  higher  ionization  energies  will  charge  transfer 
to  those  of  lower  ionization  energy.  A listing  of  the  ionization  energies  of  species 
of  interest  is  given  in  Table  2.  As  can  be  seen,  CO2  has  the  lowest  ionization 
energy  of  the  four  test  gases  of  interest  and  thus,  if  the  ion/molecule  reaction 
rates  are  sufficiently  rapid,  CO2  will  become  the  dominant  positive  ion  in  all 
these  mixtures,  (This  would  not  be  the  case  if  a significant  amount  of  O2  were 
created  while  the  beam  is  on. ) 


Table  2.  Ionization  Energies 


The  nuH-haniritiis  by  which  this  chai’sc  cxchanf'c  t)CCur.->  can  |>i'<iducc  additional 

neutral  .spocieri  and  are  discussed  below.  'I'he  rate  constants  used  in  this  discus - 

7-9 

Sion  are  taken  for  the  most  part  from  ' other  sources,  and  the  characteristic 
reaction  times  listed  are  for  the  case  of  !>0  torr,  P(’().,  0.5  torr.  (Tliese 

characteristic  times  scale  inversely  with  pressure  or  pressure  squared  for  binary 
and  tertiary  reactions  respectively. ) The  primary  reaction  mechanisms  for  each 
beam  produced  ion  are: 

(1)  N^' 

N.,'  I CO.,  - CO.,'  t N,,  k 5=  lO"^  cm'\/sec,  r i>  X lO’**  sec 

(Ki) 

N.,'  t N.,  t M > I M k « 9 X 10  cm*’/sec.  t - 7 X 10  ''  sec 

(17) 

+ CO2  * ' *^2  T <>  X 10  sec 

(18) 

C02^  • '^2  ' ‘^‘^^2  ■*  ' *^^^2  ' ^2  ~ cin^/sec,  7 ti  X 10'‘*  sec 

(19) 

(2)  n' 

n'  + CO2  * CO,,^  t N k « 1.  .'1  X 10  ^ cm'^/sec,  t •')  X 10  **  sec 

(20) 

I I -29  (!  -!! 

N + N2  ' M - N.j  ♦ M k = :l  X 10  ■ cm  /sec,  t 1.2  X 10  sec 

(21) 

N.j'  t CO2  ■■  Co.,'  ■ N2  t N k s lO"'*  cm'^/sec,  7 (>  X 10  **  sec 

(22) 

CO.,'  • N2  t Co.,  . CO,,'  • C(l,  I N.J  k - lO'*’  cm'Vsi<^f.  T H X lO'*'  sec 

(19) 


7.  Pehsenfeld,  P.C.,  .Schmeltekopf,  A.l..,  Ounkin,  I).  11. , and  b'erfiuson,  I'.  M. 

(19(>9)  Compilation  of  Reaction  Hate  Const. mts  Measured  in  the  ICS.SA  Plow- 
ing Afterglow  System  to  August  1909,  IIS.SA  Technical  lleport  lOh’l,  I'l'i-Al.O 

8.  Ferguson,  F.  F.  (1978)  Atom.  Pat,  and  ,\ucl.  I)at.  Tables  12:l.'i9. 

9.  flood.  A,  (1979)  Chem.  llev,  7.'>;r)01. 


<3)  O'^ 


+ CO,  - + CO  k = lO"®  cm^/sec,  t = 6 X 10  ® sec 

£ A 

(23) 

(4)  CO^ 


co"^  + COg  - + CO 

k = 1.  1 X lO'®  cm^/sec,  t = 6 X 10  ® sec 

(24) 

(5)  C'^ 

C^  + COg  - co"^  + CO 

k = 1.  9 X 10*^  cm^/sec,  t = 3 X 10  ® sec 

CO^  + COg  - COg^  + CO 

(25) 

k = 1.  1 X lO'®  cm^/sec,  t = 6 X 10  ® sec 

(24) 

As  can  be  seen,  with  the  exception  of  O^,  all  directly  produced  ionic  species 
charge  transfer  to  COg^.  Ultimately,  the  two  terminal  ions  COg  and  Og  recom- 
bine producing  neutral  species.  The  appropriate  reaction  mechanisms  are 

COj^  + e - CO  + O 

(26) 

02^  + e - O + O 

(27) 

and  therefore  the  ultimate  conversion  of  ionic  to  neutral  species  may  be  repre- 
sented by  the  global  reactions 


N2^  + e + COg  - N2  + CO  + O 

(28) 

+ e + COg  - N + CO  + 0 

(29) 

O"*^  + e + COg  - CO  + 20 

(30) 

00"*^  + e + COg  - 2CO  + 0 

(31) 

C"^  + e + 2CO2  - 3CO  + 0 . 

(32) 
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I 

From  consideration  of  the  global  reactions  (28)-(:i2),  the  fractional  ionization 
levels  of  Table  1,  and  the  dissociation  cross  section  of  N,  shown  in  Figure  7,  the  | ' 

total  production  of  neutrals  per  ion  pair  may  be  specified.  These  are  listed  in 
Table  3.  From  the  analysis  above,  species  such  as  C,  NO,  CN,  and  O^  are  not 
created  (except  possibly  in  trace  quantities),  through  the  sequence  of  electron  or 
ion  reactions  occurring  in  the  gas.  Furthermore,  the  dominant  source  of  CO  pro- 
duction, reaction  (2(>),  has  been  identified. 

I 

1 


H 


Similar,  although  less  complicated,  analysis  can  be  provided  for  the  Ar/C02 
and  He/C02  mixtures.  In  the  case  of  Ar,  the  dominant  ion  produced  is  Ar^  and 
the  relevant  reaction  sequence  is 

k = ()  X 10  cm'^/sec,  t 1 X 10  ^ sec 

(33) 

k 2X10'^'^  cm*’/sec,  t ' 1.8X10  *’  sec 

(34) 

k 1.  1 X 10  cm'^/sec,  T = (i  X 10  ^ sec 

(3r,) 

and  thus  each  Ar*^  produces  a C02^  which  upon  recombination  produces  one  CO 
molecule  and  one  oxygen  atom  per  ion  pair.  The  case  of  Helium  is  sornewliat 
more  complicated  because  the  reaction  between  He*  and  COg  has  three  branches 


Ar^  f CO2  ^ C02^  Ar 
Ar^  + Ar  + M -*  AK2'  + M 
Ar2^  + CO2  ♦ ^^^2^  ’ ^ 


Table  3.  Species  Production/lon  Hair 


Production  Per 

Mixture 

Species 

Ion  Pair 

N2/CO2  100 

N 

1.  (> 

0 

1.0 

CO 

1.0 

Ar/C02  - 100 

0 

1.0 

CO 

1.  0 

He/COg  - 100 

0 

2.  1 

CO 

2.  1 

i* 
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k = . 8 X 10  cm^/riec,  t ' 7 X 10  **  see 

CIO) 

k = 1.  5 X 10  cm'^/sec,  t 4 X 10  ^ .see 

(87) 

k 5 X 10  em'^/,sec,  t - 1.2X10  sec  . 

(38) 

I The  CO*^  created  will  charge  exchange  to  CO2  via  reaction  (24);  however,  the  o' 

will  form  the  terminal  ion  O.,^  via  reaction  (23).  The  net  species  production  rates 
from  this  sequence  of  reactions  is  listed  in  Table  .3.  It  is  to  be  noted  that,  unlike 
the  N.,/C02  and  Ar/CC),,  mixtures,  more  than  half  of  the  CO  created  in  Ile/CO^ 
mixtures  is  the  result  of  charge  exchange  reactions.  The  recombination  reaction 
(20)  is,  however,  the  only  mechanism  for  CO  formation  which  is  common  to  all 
three  mixtures. 

It  is  evident  that  in  all  three  mixtures,  significant  quantities  of  CO  and  free 
radicals  are  created  by  the  electron  beam  radiation.  Although  it  has  been  found 
that  electron  and  ion  reactions  occur  on  a very  short  time  scale  — and  thus  will 
only  take  place  while  the  beam  is  on  — reactions  between  neutral  species  can  occur 
on  longer  time  scales.  Furthermore,  the  concentration  of  beam-produced  neutral 
species  with  the  chamber  can  increase  from  pulse  to  pulse.  Indeed,  although  the 
e-beam  pulse  duration  is  ■ a msec  tne  gas  residence  time  within  the  test  chamber  has 
been  measured  to  be  10-15  sec  for  the  cases  shown  in  h'igures  1-4.  Thus  the  cem- 
centration  of  beam-produced  species  could  build  up  to  sufficiently  large  levels  in 
the  test  chamber  so  that  they  could  effect  the  relaxation  behavior  of  the  observed 
fluorescence. 

To  provide  an  estimate  of  this  effect  the  amount  of  energy  supplied  to  the  gas 
by  the  electron  beam,  and  the  gas  volume  to  which  this  energy  is  supplied,  must 
be  specified.  As  mentioned  earlier  the  energy  deposition  is  gas  specific,  spatially 
nonuniform,  and  will  depend  upon  the  thickness  of  the  target.  Beam  energy  degra- 
dation is  generally  defined  in  terms  of  the  electron  stopping  power,  1 /p  dK /dx, 
which  is  a function  of  beam  voltage.  The  total  energy  lost  per  electron  by  a beam 
propagating  a distance  f through  a target  of  density  p is  defined  by 

/ 

p ( i dK/dx)  dx  (30) 


He'  + CO.,  . Co'  ( O I He 
. O*^  + CO  + He 
- CO2'  ‘ lie 
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where  the  fact  that  beam  voltage  decreases  with  increasing  propagation  distance 
is  implicity  included  in  the  evaluation  of  the  integral.  Values  of  1/p  dE/dx  for  a 
beam  voltage  of  32  keV  are  listed  in  Table  4 for  the  gases  of  interest.  Also  listed 
are  gas  molecular  weights,  the  effective  range  over  which  a beam  with  initial 
energy  of  32  keV  would  be  dissipated  and  w the  average  number  of  eV  required  to 
produce  an  ion  pair  in  the  gas.  The  listed  values  of  range  and  electron  stopping 
power  were  taken  from  the  tabulation  of  Berger  and  Seltzer,  and  the  energies 
per  ion  pair  were  taken  from  the  work  of  Whyte.  The  variation  in  beam  stopping 
power  with  beam  voltage  is  also  shown  in  Figure  8 for  the  gases  Ng,  He,  and  Ar. 

Such  tabulations  are  sufficient  to  calculate  the  energy  transfer  and  ion  crea- 
tion rates  for  the  cases  of  interest.  The  integrating  sphere  used  in  the  measure- 
ments has  a diameter  of  =16  cm  and  thus  AE  may  be  determined  from  Eq.  (39) 
using  the  electron  stopping  powers  plotted  in  Figure  8.  The  power  supplied  to  the 
gas  when  the  beam  is  on,  is  then’defined  by 

P = lAE  (40) 

and  the  rate  for  creation  of  ion  pairs  within  the  test  chamber  is  defined  by 

= P/we  (41) 

where  e is  the  unit  electron  chp-ge.  Furthermore,  the  total  energy  supplied  to 
the  gas  per  pulse  and  number  of  ion  pairs  created  per  pulse  may  be  determined  by 
multiplying  Eqs.  (40)  and  (41)  by  the  pulse  duration  t^.  Finally  one  must  evaluate 
the  total  energy  transferred,  and  ion  pairs  created,  during  the  gas  residence  time 
within  the  test  chamber.  These  quantities  may  be  determined  from  Eq.  (40)  and 
(41)  by  multiplying  by  the  product  of  the  gas  residence  time  and  the  beam  duty 
cycle  n.  These  quantities  are  listed  in  Table  5 for  the  three  mixtures  of  interest. 
The  parameters  used  in  evaluating  these  quantities  were  I = 0.  6 mA,  4.  6 msec, 

Tp  ^ 10  sec  and  rj  - 0.  2. 

The  quantities  listed  in  Table  5 are  essentially  integrals  of  the  total  energy 
and  number  of  ion  pairs  created  within  the  two  liter  integrating  sphere.  It  would 
be  valuable  to  determine  the  volume  of  gas,  within  the  sphere,  which  is  irradiated 
by  the  beam  so  that  the  local  densities  of  beam  produced  species  could  bo  defined. 
This  determination  is  complicated  by  the  fact  that  the  beam  energy  deposition  is 
spatially  nonuniform.  For  example  the  Ar  and  N2  mixtures  may  be  considered  as 
optically  thick;  while  the  electron  range,  as  tabulated  in  Table  4,  being  28  cm.  A 

10.  Berger,  M.J.,  and  Seltzer,  S.M.  (1964)  Tables  of  Energy  Bosses  and  llanges 

of  Electrons  and  Positions,  NASA  SP-3012. 

11.  Whyte,  G.  N,  (1963)  Radiation  lies.  18;2r)r). 
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Table  4.  Electron  Degradation  Propertie; 


Molecular 


Stopping 

Power 

MeV  -cm2/gm 


Energy  l^er 
Ion  Pair 


gm/mole 


I'able  5.  lUectron  IJoam  Gas  Interaction  Parameters 


Gas  Mixture 

50  torr  lle- 
0,  5 torr  CC).^ 

- .. 

50  torr  N2- 
0.  5 torr  CO2 

50  torr  Ar- 
0,  5 torr  CO.^ 

Density, 
gms  'cc 

1.2  X 10‘'’ 

7.  5 X 10"^ 

1. 08  X 10"'* 

AE,  eV 

1.  8 X 10^ 

1.  1 X lO'^ 

0 

X 

P watts 

1.  1 

G.  G 

G.  G 

' -1 

N . sec 

17 

1.  fi  X 10 

I. 2 X 10^^ 

CO 

0 

X 

Energy/ 

Pulse 

Joules 

5,  1 X lO"'^ 

3 X 10'^ 

3 X 10"^ 

Ion  Pairs/ 

Pulse 

14 

7.4  X 10 

5.  5 X lo’-^ 

1 5 

7.4  X 10 

Energy  in  t„ 
Joules 

2,  2 

13.  2 

13.2 

Ion  Pairs  in 

■^R 

17 

3.  2 X 10 

2.4  X 10^*^ 

1 Q 

3.  2 X 10 

measurement  of  the  volume  energy  distribution  rate  in  a thick  target  is  shown  in 

19 

Figure  9,  as  excerpted  from  Cohn  and  Caledonia.  " As  can  be  seen  both  the  axial 
and  radial  variation  of  the  deposition  rate  is  quite  significant  over  the  cloud  length. 
Indeed  it  can  be  seen  that  variations  of  over  an  order  of  magnitude  in  local  ion 
pair  creation  rates  and  energy  deposition,  will  occur  within  the  test  chamber. 
Although  the  helium  mix  may  be  considered  a thin  target,  the  energy  deposition  in 
this  case  will  also  vary  significantly  over  the  length  of  the  integrating  sphere  (thin 
target  behavior  will  be  discussed  in  more  detail  later  in  the  text). 

Thus,  considerable  effort  would  be  required  in  order  to  define  the  pulse-to- 
pulse  beam  energy  deposition.  An  approximate  estimate  can  be  provided  for  the 
long  time  build  up  of  beam-produced  species  and  energy  because  of  diffusive  effects. 

A typical  value  of  the  diffusion  coefficient  for  species  such  as  CO  diffusing  in  or 

2 

Ar  at  pressures  of  50  torr  would  be  ~3  cm  /sec.  Atomic  species,  or  species  in  a 
helium  carrier  gas  would  exhibit  even  larger  diffusion  coefficients.  The  charac- 
teristic distance  L in  which  species  will  diffuse  in  a time  t may  be  approximated 
by  the  relationship 


12.  Cohn,  A,,  and  Caledonia,  G.  (1970)  J.  Appl.  Phys.  4 1 :37(i7. 
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Figure  9.  Normalized  Contour  Plot  of  Volume  Energy 
Deposition  Rate  in  a Thick  N,  Gas  Target.  L is  the  elec- 
tron practical  range,  (excerpted  from  ref.  12) 

L~  (Dt)^/^  (42) 

where  D is  the  diffusion  coefficient.  For  the  gas  residence  time  of  10  sec,  this 
characteristic  dimension  would  be  *5.  5 cm,  which  is  comparable  to  the  8 cm  rad- 
ius of  the  integrating  sphere.  Thus,  to  first  order  the  volume  of  beam  — effected 

3 3 

gas  may  be  estimated  to  be  10  cm  . From  Table  5,  that  would  imply  that  the 

3 

average  number  of  ion  pairs  produced  per  cm  during  the  gas  residence  time  would 
be  *3.2X10^^,  2.4X10^^,  and  3.  2 X 10^^  for  the  He,  N2,  and  Ar  mixtures, 
respectively.  Furthermore,  the  energy  imparted  to  the  gas  by  the  beam  would 
cause  an  increase  in  gas  temperature.  The  estimated  gas  temperature  increase 
over  the  residence  time  would  be  40,  160  and  240°K  for  the  He,  N2,  and  Ar  mixes, 
respectively. 

It  should  be  remembered  that  these  estimates  are  crude  and  could  be  off  by  a 
. factor  of  two  or  more  because  of  uncertainties  in  the  definition  of  the  volume  of  gas 

affected  by  the  beam  and  in  the  gas  residence  time.  Nonetheless,  these  perturba- 
tions are  large  and  most  probably  would  play  a role  in  the  kinetics  controlling  the 
observed  vibrational  fluorescence.  Specifically  from  Table  3,  it  can  be  seen  that 
for  each  ion  pair  formed  at  least  one  oxygen  atom  and  CO  molecule  are  produced 
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and,  in  the  case  of  nitrogen,  1.  6 nitrogen  atoms  are  also  produced.  Since  it  has 
been  shown  that  3 X 10^^  - 3 X 10^^  ion  pairs  are  produced  per  cm^  during  the  gas 
residence  time,  then  at  least  that  many  O,  CO  and,  depending  on  mixture,  N will 
ultimately  be  produced.  These  concentrations  are  to  be  compared  to  the  initial 
CO2  concentration  of  1.  6 X 10^^/cm^. 

The  actual  concentration  of  beam  produced  neutral  species  need  not  increase 
to  such  high  concentrations  inasmuch  as  neutral  species  chemical  reactions  will 
occur  within  the  test  chamber.  These  chemical  reactions  could  lead  to  the  recom- 
bination of  beam  produced  neutral  species  or  alternately  to  the  creation  of  new 
species.  Thus,  the  neutral  chemistry  occurring  within  the  chamber  must  be 
modeled  in  order  to  determine  the  average  composition  of  the  electron  irradiated 
gas  within  the  chamber.  The  dominant  neutral  kinetic  mechanisms  occurring  in 
a mixture  of  No/CO„  dosed  with  O,  N,  and  CO  are  listed  in  Table  G.  The  rate 

i c 13  . 

constants  for  these  reactions  were  taken  from  Garvin  and  Hampson  with  the 

14 

exceptions  of  reactions  (d)  and  (h)  which  were  deduced  from  the  studies  of  Inn 

15  1 G 

and  Slanger  et  al,  ^ and  from  Rawlins  and  Kaufman,  respectively. 

Reaction  (h)  of  Table  6 is  perhaps  the  most  crucial  reaction  in  the  sequence. 

For  example,  if  the  reaction  between  nitrogen  atoms  and  CO2  were  rapid,  then  an 

additional  CO  molecule  could  be  produced  for  every  nitrogen  atom  created  while 

the  electron  beam  is  on.  Furthermore  alternate  reaction  paths,  which  would 

ultimately  create  species  such  as  CN,  might  be  possible.  Fortunately,  while  this 

1 6 

analysis  was  in  progress,  Rawlins  and  Kaufmann  published  a timely  study  on  the 

reaction  of  CO2  with  active  nitrogen.  In  this  work,  it  was  found  that  the  dominant 

reaction  between  nitrogen  atoms  and  CO,  occurred  in  a several  step  mechanism 

^ * 

involving  an  excited  electronic  state  of  nitrogen,  N2  • The  general  reaction 
sequence  proposed  by  Rawlins  and  Kaufmann  to  describe  their  observations  is 

kM 

N + N + M 5 N2*  + M 
k’^ 

N2*  + M -9  N2  + M (43) 


13.  Garvin,  D. , and  Hampson,  R.F. , editors  (1974)  Chemical  kinetics  data 

survey  VII,  Tables  of  Rate  and  Photochemical  Data  for  Modeling  at  the 
Stratosphere  (Revised),  NBSIR  74-430. 

14.  Inn,  E.C.Y.  (1974)  J.  Chem.  Phys.  61:1589. 

15.  Slanger,  T.G.,  Wood,  B.J.,  and  Black,  G.  (1972)  J,  Chem.  Phys.  57:233. 
IG.  Rawlins,  W.T,,  and  Kaufman,  P.  (1976)  J.  Chem.  Phys.  64:1 128. 
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3 6 

k {crn  /sec  or  cm  /sec) 

(a) 

O + O + Ng  - Og  + Ng 

4.8  X 10'^^  (300/T) 

(b) 

O + N + Ng  ->  NO  + Ng 

1.  8 X 10"^^ 

(c) 

O + Og  + Ng  - O3  + Ng 

1.05X10-34 

(d) 

0 + CO  + Ng  - COg  + Ng 

0.8X10-33 

(e) 

0 + O3  - 2 Og 

1.9X10-11  e-3300/T 

(f) 

N + N + M-.  Ng  + M 

8.3X10-34  e300/T 

(g) 

N + NO  ^ Ng  + 0 

2.7  X 10-11 

(h) 

N + N + COg  - CO  + Ng  + 0 

N + COg  - NO  + CO 

jsee  text 

(i) 

O3  + Ng  ^ 0 + Og  + Ng 

6.6X10-1°  e-ll-4°°/T 

(j) 

N + Og  - NO  + 0 

l.lX10-14Te-315°/T 

(k) 

N + O3  - NO  + Og 

5.  7 X 10-13 

(1) 

0 + NO  + Ng  - NOg  + Ng 

4.  2 X 10-33 

(m) 

0 + NOg  - NO  + Og 

9.  1 X 10-13 

(n) 

0 + NOg  + Ng  NO3  + Ng 

1.  X 10-31 

(o) 

O3  + NO  - NOg  + Og 

9 X 10-13  e-13°°/T 

(P) 

O3  + NOg  ^ NO3  + Og 

1.  IX 10-13  e-3450/T 

(q) 

N + NOg  ^ 2 NO  ( ?) 

1.  85  X 10-11 

(r) 

NO  + NO3  ^ 2 NOg 

^10-11 

I 


This  rate  is  based  on  the  authors  interpretation  of  studies  presented  in  refs.  14 
and  15. 


■where,  for  the  present  case,  M represents  N2  and  CO^-  (This  sequence  of  reac- 
tions is  further  complicated  in  that  Ng  may  be  produced  in  several  different  vibra- 
tional states  each  of  which  has  a different  quenching  constant. ) The  reaction 
sequence  (43)  may  be  represented  by  the  global  reaction 
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N + N + COg  " N2  + CO  + O 


(44) 


where  the  rate  constant  is  a function  of  the  ratio  of  N2  to  CO 


It  may  be 

determined  from  the  work  of  Whyte,  that  a value  k^^  = 4 X 10 cm  /sec  is 
appropriate  for  the  conditions  of  the  present  analysis.  Thus  the  only  important 
neutral  reaction  occurring  in  the  test  chamber  which  involves  CO2  a relatively 
slow  three  body  process.  From  examination  of  Table  G it  can  be  seen  that  the 
same  is  true  regarding  reactibns  involving  CO. 

E.  Lee  of  AFGL  has  exercised  the  reaction  scheme  of  Table  G in  order  to 
determine  the  final  neutral  species  products  produced  in  the  test  chamber.  His 
calculation  included  pulsed  source  terms  for  O,  CO,  and  N as  per  the  beam  duty 
cycle.  It  was  found  that  on  time  scales  less  than  1 sec  the  oxygen  atom  concentra- 
tion approaches  "steady  state"  and  the  two  dominant  species  increasing  in  time 
without  significant  loss  are  CO  and  O2.  At  larger  times,  species  such  as  NO  and 
Og  begin  to  play  a role  in  the  kinetics.  It  will  be  recalled  that  O2  was  not  consid- 
ered as  a collision  partner  in  the  previously  discussed  analysis  of  the  ionic  kine- 
tics. If  sufficient  O2  is  present  in  the  chamber,  O2  w’ill  become  the  dominant 
ion  and  the  production  of  vibrationally  excited  CO  via  electron  ion  recombination. 


reaction  (26),  would  cease.  The  primary  initial  ions  created  by  the  beam,  such  as 
N2^  and  will  not  preferentially  form  02^  over  C02^  unless  the  O2  concentration 


^2' 


CO2  are  typically  larger  than  these  for  charge  exchange  to  O,.  However,  the  CO2 
itself  will  directly  charge  exchange  to  O2  via  the  reaction 


COt  + Or. 


0„  + COr 


k = 5 X 10  cm'^/sec 


(45) 


Thus  when  the  O2  concentration  is  sufficiently  large  so  that  reaction  (45)  will  occur 


more  rapidly  than  electron-ion  recombination,  the  dominant  positive  ion  will  be 
O,^  rather  than  CO„^.  The  characteristic  time  for  recombination  is  estimated  t 
be  10  psec  and  thus  the  condition  for  O2  to  become  the  dominant  ion  becomes 


5 X lO"*’  N 


O, 


> 10' 


(4  6) 


;)o 


F’ 


or 


(47) 


H 


i»2X10^^/cm^ 

It  can  be  seen  from  reference  to  Table  5 that  the  density  of  oxygen  molecules 
listed  above  is  of  the  same  order  of  magnitude  as  that  predicted  to  build  up  within 
the  chamber  during  the  gas  residence  time.  Thus  the  production  mechanism  for 
vibrationally  excited  CO  will  be  partially  effected  by  the  presence  of  beam-produced 
species. 

The  arguments  presented  here  for  N2/CO2  mixtures  will  also  be  true  for  the 
Ar/C02  mixture.  The  neutral  reaction  mechanism  in  the  latter  mixtures  will  be 
the  same  as  that  listed  in  Table  6 with  the  exception  of  the  reactions  involving  N 
and  N2.  In  the  case  of  the  He/C02  mixture,  CO  and  O2  will  again  be  the  dominant 
beam  produced  neutrals;  however,  the  0„  concentration  will  not  be  sufficiently  high 

f ‘ 

to  effect  the  CO2  concentration. 

The  conclusion  of  the  kinetic  analysis  is  that,  over  the  gas  residence  time  of 
10  sec,  significant  quantities  of  electron  beam  produced  neutral  species  will  build 
up  in  the  test  chamber  for  the  gas  and  electron  beam  conditions  typical  of  the  data 
base.  These  beam  produced  species  will  be  dominantly  O2  and  CO  with  traces  of 
O,  Og  and  in  the  case  of  the  nitrogen  mix,  NO  and  N02*  Indeed  the  predicted  con- 
centrations of  CO  and  O2  for  the  N2  and  Ar  mixes  are  within  10-20  percent  of  the 
initial  CO2  concentrations  (this  of  course  implies  that  the  CO2  concentration  has 
decreased  by  this  much).  Furthermore,  significant  increases  in  gas  temperature 
have  also  been  predicted.  These  phenomena  effectively  preclude  the  determination 
of  any  kinetic  data  from  the  observed  CO  fluorescence  histories;  however,  the 
analysis  developed  may  be  used  to  define  the  experimental  conditions  required  in 
order  to  take  viable  data. 

Before  proceeding  with  this  definition,  it  is  worthwhile  reviewing  the  available 
rate  constant  data  base  for  the  kinetic  mechanisms  of  interest  in  order  to  ascertain 
the  maximum  acceptable  levels  of  beam  produced  species  within  the  test  chamber. 
The  basic  kinetic  mechanisms  to  be  studied  are 

CO(v)  + N2(v  - 0)  - CO(v  - 1)  + N2(1)  (2) 

CO(v)  + CO2  - CO(v  - 1)  + C02(v3)  . (3) 

The  available  room  temperature  rate  constants  for  these  two  slightly  endothermic 
processes  are  shown  in  Figure  10.  The  rate  constants  for  processes  (2)  and  (3) 
for  V = 1 are  well  defined.  The  value  shown  in  Figure  10  for  N2  is  representative 
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CO  Vit)rationdl  level 


Figure  10.  Room  Temperature  Rate  Constants 
for  the  Reactions  CO(v)  + CO2  -*  CO(v  - 1)  + 
COgCvg)  and  CO(v)  + Ng  - CO(v  - 1)  + N2(v  = 1) 


17-19 

of  measurements  provided  in  references  and  the  value  shown  for  CO,  is  as 

20-23  ^ 

reported  in  references.  The  data  shown  for  vibrational  levels  ^ for  both 

24 

N2  and  CO2  are  taken  from  the  work  of  Hancock  and  Smith.  No  measurements 
are  available  for  v = 2,  3.  Note  the  decreasing  rate  constant  with  increasing 
vibrational  level  is  an  expected  behavior  since  the  anharmonicity  of  the  CO  mole- 
cule is  such  that  the  endothermicity  of  reactions  (2)  and  (3)  increases  with  increas- 
ing vibrational  level.  The  observed  increasing  rate  constant  process  (3)  at 


17.  Green,  W.  H. , and  Hancock,  J,  K.  (1973)  IEEE  J.  Quantum  Electr.  QE-9;50. 

18.  Starr,  D.  F. , Hancock,  J.K.,  and  Green,  W.  H.  (1974)  J.  Chem.  Phys.  61: 

5421.  

19.  Zittel,  P.F.,  and  Moore,  C.B.  (1972)  Appl.  Phys.  Letters  21:81. 

20.  Rosser,  W.A.,  Jr.,  Sharma,  R.D.,  and  Gerry,  E.T.  (1971)  J.  Chem.  Phys. 

M:1196. 

21.  Stephenson,  J.C.,  and  Moore,  C.B.  (1972)  J.  Chem.  Phys.  56:1295. 

22.  Miller.  D.J,,  and  Millikan,  R.C.  (1974)  Chem.  Phys.  6:317. 

23.  Starr,  D.  F. , and  Hancock,  J.  K.  (1975)  J.  Chem.  Phys.  63:4730. 

24.  Hancock,  G. , and  Smith,  I.  W.  M.  (1971)  Appl.  Opt.  10:1827. 
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vibrational  levels  2 7 is  surprising  and  it  has  been  postulated  to  occur  because 
of  a contribution  from  the  reaction 

CO(v)  + CO2(01^0)  - CO(v  - 1)  + CO2(0lh)  . (48) 

The  other  mechanism  involving  CO  and  the  test  gases  is  vibration  to  translation 
(V-T)  deactivation 

CO(v)  + M - CO(v  - 1)  + M (49) 

where  M can  be  Ng,  CO2,  Ar  or  He.  The  vibrational  level  dependent  rate  constants 

for  process  (49)  have  been  measured  only  for  helium  and  are  shown  in  Figure  11. 

25 

The  rate  constant  for  v = 1 is  from  the  work  of  Miller  and  Millikan  while  those 

24 

for  V - 9 - 13  are  from  Hancock  and  Smith.  Also  shown  is  a theoretical  predic- 
tion for  the  rate  constant  for  process  (49)  with  M = N„.  This  prediction,  taken 

26  27  ^ 

from,  is  based  upon  SSH  Theory  and  is  representative  of  the  expected  behavior 

and  magnitude  of  the  V-T  rates  for  M = CO2,  Ar.  Note  — from  a comparison  of 
the  data  in  Figures  10  and  11  — it  can  be  seen  that  CO  V-T  deactivation  by  He  will 
be  competitive  with  CO  VV  exchange  to  CO2  in  100:1  mixtures  of  He-C02.  On  the 
other  hand,  V-T  collisions  with  N2,  CO2  and  Ar  will  be  unimportant. 

The  dominant  beam  produced  neutral  species  will  be  vibrationally  cold  CO 
which  can  react  with  vibrationally  excited  CO  via  the  intramode  vibration  to  vibra- 
tion (V-V)  exchange  reaction 

CO(v)  + CO(0)  - CO(v  - 1)  t cod)  . (50) 

Process  (50)  does  not  decrease  the  number  of  vibrational  quanta,  but  rather  rear- 
ranges the  distribution  of  quanta,  and  thus  can  effect  the  relaxation  of  individual 

24  2fi  '^9 

vibrational  levels.  Several  ’ ’ measurements  of  the  room  temperature  rate 

constants  for  process  (50)  are  shown  in  Figure  12  and  they  can  be  seen  to  be  in 
agreement  to  within  a factor  of  two.  The  rate  constant  for  v 1 is  not  shown  since 
it  corresponds  to  a null  reaction.  The  maximum  concentration  of  cold  CC  which 
can  be  allowed  in  the  system  may  be  deduced  by  comparing  the  rate  constants  of 

25.  Miller,  D.J.,  and  Millikan,  Fl.C.  (1970)  J.  Chem.  Phys.  M:3384. 

28.  Center,  R.E.,  and  Caledonia,  G.  F.  (li)71)  Appl.  Opt.  10:1795. 

27.  Schwartz,  R.N.,  Slaw,-,ky,  Z. , and  Horzfeld,  K.  F.  (1952)  J.  Chem.  Phv-.. 

20:1591.  ■ 

28.  Liu.  y.  S.  . McFarlane.  R.A.,  and  Volga,  G.J.  (1975)  J.  Chem.  Pliys. 

W:228. 

29.  Fushiki,  Y. , and  Tsuchiya,  S.  (1974)  Jap.  J.  Appl.  Phys.  13:1043. 


Figure  11.  Rate  Constants  for  Vibration  to  Trans- 
lation Deactivation  of  CO  ( O 24,  —ref.  2R) 
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Figure  12.  Room  Temperature  Rate  Constants  for 
the  Reaction  C'Oiv)  f CO(0)  • COfv  - 1)  ‘ CO(l) 

(O  fef.  24,  O ref.  2H,  0 ref.  2f)) 


process  (50),  Figure  12,  with  those  for  reactions  (2)  and  (3),  Figure  10.  For 
example,  if  it  is  required  that  reaction  (50)  have  no  more  than  a 10  percent  effect 
on  the  relaxation  of  any  vibrational  level  then  it  is  required  that  [CO]  < 10  [CO2I 
or  10”^  [Njl . For  a 1 percent  effect,  the  CO  concentration  would  have  to  be  an 
order  of  magnitude  smaller  than  those  listed  above.  As  a numerical  example,  for 
the  Ar/C02  mix  considered  earlier  the  CO2  concentration  was  1.  6 X 10  /cm  and 
thus  - if  it  is  desired  that  cold  CO  contribute  less  than  10  percent  to  the  observed 
relaxation  of  any  vibrational  level  — the  CO  concentration  must  be  less  than 
1.  6 X 10^^/cm^.  This  CO  density  is  two  orders  of  magnitude  lower  than  the  pre- 
dicted density  of  the  preliminary  measurements. 

The  only  other  diatomic  molecule  created  in  large  quantities  by  the  electron 
beam  is  02*  This  species  may  react  with  CO  via  the  V-V  exchange  reaction 

CO(v)  + O2  - CO(v  - 1)  + 02(v  = 1)  . (51) 

24  30 

The  room  temperature  rate  constant  for  process  (51),  v = 1,  has  been  found  ’ 

-17  3 

to  have  the  low  value  of  6 X 10  cm  /sec.  However,  since  the  reaction  is  exo- 
thermic, the  rate  constant  is  expected  to  increase  with  increasing  CO  vibrational 
level;  that  is,  the  exothermicity  decreases  with  increasing  v and  thus  the  reaction 
approaches  energy  resonance.  Indeed  Hancock  and  Smith  have  determined  the 
rate  constants  for  process  (51),  v = 12,  13,  to  be  3.  1 X 10  ^ and  4.  4 X 10 
cm^/sec,  respectively.  However,  these  rate  constants  are  still  lower  than  the 
respective  values  for  CO2.  and  thus,  since  the  concentration  of  beam  produced  O2 
must  be  less  than  one  half  that  of  CO,  the  previous  limits  specified  for  the  CO 
concentration  will  be  more  than  ample  in  precluding  any  undesired  relaxation  due 
to  reaction  (51). 

The  last  consideration  is  to  the  possible  V-T  deactivation  of  CO  by  oxygen  and 
nitrogen  atoms.  Little  kinetic  information  is  available  on  these  processes. 

O 1 

Center  has  measured  the  rate  constant  for  V-T  deactivation  of  CO  (v  1)  by 
O-atoms  at  high  temperatures.  An  extrapolation  of  his  I'csults  to  room  tempcratur 
would  result  in  a rate  constant  of  10  ' cm  /sec.  It  would  be  expected  that  tliis 

rate  constant  would  increase  with  increasing  vibrational  level,  perhaps  being  a 
few  orders  of  magnitude  larger  by  v 13.  Since  the  oxygen  atom  concentration 
will  perforce  be  less  than  or  equal  to  the  CO  concentration,  the  previously  defined 
CO  density  limit  would  also  seem  to  he  sufficient  to  rule  out  any  significant  vibra- 
tional relaxation  resulting  from  the  presence  of  oxygen  atoms. 


30.  Miller,  D.J.,  and  Millikan,  It.C.  ( 1974)  ClieTn.  I'hy-..  l.ettcr-.  27:10. 
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In  conclusion,  it  was  found  that  the  experimental  conditions  have  to  be  recon- 
figured to  reduce  the  concentration  of  beam  produced  species  within  the  test  cham- 
ber by  at  least  two  orders  of  magnitude.  This  could  be  achieved  in  a number  of 
ways  (a)  decrease  the  gas  residence  time,  (b)  decrease  the  beam  duty  cycle, 

(c)  work  with  thin  rather  than  thick  gas  targets,  and  (d)  decrease  the  beam  current. 
Several  of  these  steps  were  taken  simultaneously  in  order  to  develop  a new  data 
base  which  would  provide  the  desired  kinetic  information.  The  analysis  of  this 
data  is  described  in  the  next  two  sections. 


:i.  M)  ii\r\ 

A number  of  steps  were  taken  in  order  to  attain  the  desired  experimental 
operating  conditions.  The  most’  prominent  of  these  were  the  remov^al  of  the  gold 
coated  integrating  sphere,  which  had  provided  the  limiting  impedance  to  the  gas 
flow,  and  the  opening  of  a secondary  pumping  port.  A flow  meter  was  also  added 
to  the  system,  in  order  to  better  determine  the  gas  residence  time,  and  alcohol/ 
dry  ice  traps  were  implemented  to  reduce  test  gas  impurities.  A series  of  runs 
were  then  performed  in  "thin  target”  mixtures  of  Ar/C02,  this  mixture  being 
considered  the  most  straightforward  to  analyze  based  upon  the  considerations  dis- 
cussed in  -Section  2.  Fluorescence  measurements,  with  optical  resolution  of 
10  cm  \ were  performed  for  an  argon  pressure  of  9 torr  and  nominal  COg  pres- 
sures of  0.  1,  0.  2."),  0.  5,  and  1. 0 torr  respectively.  The  beam  conditions  were: 
current  of  0.  .1  to  0.9  mA,  voltage  of  30  kV,  and  pulse  time  of  2.4  msec  with  a 
duty  cycle  of  =10  percent.  The  gas  pumping  rate  within  the  test  chamber  was 
measured  to  be  five  standard  liters /min. 


3.1  Irrailialcil  (ins  l’r<i|irrlics 

The  distribution  of  beam-produced  species  for  one  of  these  cases,  1 '4 

torr  and  I 0.  9 mA,  will  be  worked  out  in  detail  in  order  to  demonstrate  that  local 

concentrations  of  cold  CO  are  sufficiently  low  so  as  not  to  interfere  with  the  Ci) 

vibrational  relaxation. 
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H.  Center  has  studied  ttie  growth  behavior  of  a well  collimated  electri'n 

beam  as  a function  of  target  thickness.  Ilis  observation,  for  the  case  of  thin  N,, 

and  Ar  targets,  was  that  the  beam  current  distribution  was  apparently  Gaussian 
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about  the  bean:  axis.  Center  defined  a parameter  o-'’ 

measured  from  beam  origin,  at  which  one  half  of  the  total  beam  cui’rent  is  encom- 
passed. lie  found  that  this  quantity  could  be  adc()uatcly  represented  by  a relation- 
ship of  the  form 

32.  Center,  H,  K.  (1970)  Phys.  Fluids  13:79. 


(52) 


1/2  ■ 1 + B(Ni) 


where  N is  gas  number  density,  t is  the  distance  from  beam  origin  and  the  quan- 
tities A and  B are  both  gas  and  beam  voltage  specific. 

The  values  of  the  quantities  A and  B,  as  determined  by  Center,  are  shown  in 
Table  7.  These  values  have  been  used  to  predict  the  expected  half-beam  radius  at 
the  end  of  the  test  chamber,  which  has  a length  of  Ki  cm.  These  predictions  are 
shown  in  Figure  13,  where  an  initial  beam  radius  of  0.  05  cm  has  been  applied. 

The  vertical  dashes  in  Figure  13  represent  the  largest  line  densities  for  which 
Center  made  measurements.  As  can  be  seen  the  beam  radius  is  largest  at  low 
voltages  and  high  line  densities,  Nf,  and  the  beam  grows  more  rapidly  in  Ar  than 
in  N2.  (Note  that  the  radii  shown  can  be  scaled  to  other  positions  in  the  test  cham- 
ber by  dividing  by  16  cm  and  multiplying  by  the  appropriate  length.)  The  total 
volume  within  the  test  chamber  irradiated  in  a single  pulse  has  been  estimated,  to 
first  order,  by  assuming  the  beam  growth  to  be  linear  (this  is  an  overestimate). 

The  resulting  volume  predictions  for  the  cases  of  Figure  13  are  shown  in  Figure  14. 
(It  is  to  be  remembered  that  this  is  only  the  volume  in  which  one  half  the  beam  cur- 
rent is  located.)  As  can  be  seen  at  low  voltages  and  moderate  line  densities,  these 

3 

volumes  can  be  several  hundred  cm  . In  comparison,  the  thick  target  "diffusion" 

3 

volume  estimated  in  Section  2.2  was  only  one  thousand  cm  . One  important 
distinction  between  the  two  cases  is,  however,  that  in  the  case  of  a thin  target, 
only  a small  portion  of  the  beam  energy  is  deposited  in  the  irradiated  volume. 


Table  7.  Beam  Growth  Parameters 


c ^ ^ — : -n 

Beam  Voltage 

A X 10^*^ 

B X 

Gas 

kV 

cm'^ 

ci‘n“ 

^2 

20 

4.  00 

2. 94 

^2 

40 

1. 00 

1.43 

■^2 

100 

0.  262 

0.  67 

Ar 

40 

2.  36 

2.  08 

A r 

100 

0.  3« 

1.  19 

37 


Bean  Half  -Current  RacJiuS  at  L ■ 16  cr^,  cn 
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Figure  13.  Electron  Beam  Half  Current  Radius 
at  L = 16  cm,  vs  Gas  Target  Line  Density 


P’ri'ssiire.  lorr 


Figure  14.  Electron  Beam  Irradiated  Volume 
in  a 16  cm  Chamber  vs  Gas  Target  Line  Density 


•t 
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An  extrapolation  of  Center's  measurements  may  be  used  to  estimate  the  beam 
dimensions  for  the  case  of  interest,  corresponding  to  a beam  voltage  of  30  kV  and 
a dominantly  argon  mix.  The  appropriate  extrapolated  constants  for  a beam  volt- 
age of  30  kV  are  A 4 X lO’"^  cm^  and  B 2.  7 X 10”^^  cm^.  The  interferometer 
is  situated  approximately  30  cm  from  the  test  chamlier  window  and  has  a circular 
field  of  vdew  of  radius  4.  1 cm  at  chamber  center.  The  center  of  the  field  of  view 
is  situated  !i  cm  downstream  of  the  electron  beam  entrance  and  thus  will  encom- 
pass the  beam  over  a distance  from  4-12  cm  downstream  of  the  beam  origin.  The 

half  beam  radius,  r,  as  defined  by  Eq.  (52)  and  a target  number  density  of 
17  3 ^ ^ 

3 X 10  /cm  , will  vary  significantly  across  the  field  of  view,  having  values  of 
0.5,  1.6,  and  3.06  cm  at  distances  of  4,  8,  and  12  cm  respectively.  This  beam 
growth  is  shown  schematically  in  Figure  15  where  the  quantities  ^^^2  ^^1/2 

are  displayed  from  beam  origin  to  the  end  of  the  field  of  view.  It  can  be  shown 
from  Figure  15  that  approximately  90  percent  of  the  beam  current  distribution 
between  3.9  and  12.2  cm  will  fall  within  the  field  of  view. 


Figure  15.  Schematic  of  Beam  Current  Distribution  Within  Field 
of  View,  Scale  1:1 
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Now  that  the  beam  current  distribution  is  defined  spatially,  local  beam  species 
production  rates  may  be  defined.  These  will,  of  course,  vary  across  the  field  of 
view  and  numerical  examples  will  only  be  provided  for  f 8 cm,  corresponding  to 
the  center  of  the  field  of  view.  The  local  ion  pair  creation  rate  per  unit  volume  a 
is  defined  by 

. i p (1/p  dE/dx)  (53) 

e w 

3 

where  p is  the  gas  density  of  2 X 10  gm/cm  , 1 Ip  dE/dx  is  the  beam  stopping 
power,  (i.  5 MeV  - cm  /gm  for  30  kV  electrons  in  Ar  (see  Figure  8),  w is  the 
average  energy  required  to  produce  an  ion  pair,  20.  1 eV  for  argon,  e is  the  unit 
electron  charge,  and  i is  the  local  current  density.  The  local  current  density  will 
be  highest  near  the  beam  axis  and  may  be  defined  there  bj' 


1/2 


Ti  r 


1/2 


(54) 


where  I is  the  total  beam  current,  which  is  0.9  mA.  Combining  Eqs.  (53)  and 
(54)  and  evaluating  results  in 


0.45  X 10 

2 

^1/2 


16 


-3  -1 

cm  - sec 


(55) 


or  o = 1.  70  X 10^^  cm  ^ - sec  ^ a.t  t =8  cm. 

The  steady  state  electron  density  is  determined  from  the  competition  betw’een 
the  production  rate  a and  the  electron  recombination  reaction 


e t 002^^  CO  + O 
and  may  be  defined  by 


(20) 


1/2 


40 


if 


I 


i 


> ■ 


3334  “7  3 

The  fate  constant  k has  been  measured'  ’ to  be  4 X 10  cm  /sec  at  room 

-1/2 

temperature  and  niay  be  expected  to  scale  with  electron  temperature  as 

Thus,  to  define  the  electron  density  one  must  estimate  the  average  steady  state 
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energy  of  secondary  electrons  produced  by  the  beam.  It  has  been  shown  that  the 
average  steady  state  energy  of  electrons  in  N2  is  = 0,  1 eV.  Similar  modeling  ^ 
for  Ar  at  a pressure  of  1 atmosphere  has  shown  that  the  average  secondary  elec- 
tron energy  is  = 3 eV.  This  high  energy  is  undoubtedly  a result  of  the  fact  that 
there  are  no  low  energy,  <10  eV,  inelastic  electron-neutral  processes  in  Ar.  It 
is  anticipated  that  the  presence  of  CO,,  in  the  Ar  mixture  will  drive  the  average 
electron  energy  to  0.  1 eV,  as  predicted  for  N.,.  Thus,  tlie  appropriate  electron- 
ion  recombination  rate  constant  to  be  used  in  Eq.  (56)  is  estimated  to  be  2.  5X10  ^ 
cm^/sec.  The  average  electron  density  at  the  center  of  the  field  of  view  is  then 

n 0,85  V 10^^ /cm'^  (57) 

e ' 

and  the  characteristic  time  for  electron  decay  upon  beam  termination  is 

T (k  n ) ^ 48  iisec  . (58) 

r r e ^ 

'I'his  characteristic  time  is  much  larger  than  the  pressure-sealed  chai-acteristic 
times  for  the  ionic  reactions  discussed  in  Section  2,  thus  ensui'ing  that  Ct is 
indeed  the  dominant  positive  ion. 

Lastly,  the  total  number  of  CO  molecules  'cm^  created  at  the  center  of  tlie 
field  of  view  during  a beam  pulse  is  defined  as 


ICO] 


max 


(59) 


1 3 

wliich  for  a pulse  time  of  2.4  msec  has  the  value  of  4.3  ' 10  “ cm  . This  densiiv 
corresponds  to  the  largest  |)c)ssil>le  concent  rat  ion  of  \dbrat  ionnlly  excited  CC  mole- 
ciilos  (at  f ti  cm)  at  any  one  time,  (liquation  (59)  reprosent.s  an  njiper  bound  on 
the  concentration  of  viluat ionally  excited  CO  both  because  reaction  (26)  need  not  be 
100  [icrcent  efficient  in  producing  excited  CO  and  also  because  smiie  degree  of  Cct 
vibrational  relaxation  will  occur  while  tiie  Ije.am  is  on.  ) 1 hi.s  jiiaxijnnm  conci  nlra- 
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CO(v)  + CO(v')  i CO(v  - 1)  + CO(v'  + 1)  (4) 

may  be  neglected.  These  reactions  can  have  large  rate  constants  when  both  v and 

26  "103 

v'  are  large.  For  example,  it  is  estimated  that  « 1 X 10  cm  /sec  for 

V = v'  = 10.  However,  as  will  be  shown,  less  than  one  twentieth  of  the  excited  CO 
molecules  reside  in  level  ten  and  thus  the  characteristic  time  for  reaction  (4) 

V = v'  = 10  is  > 50  msec.  This  is  to  be  compared  to  the  radiative  decay  time  for 
this  level  which  is  4 msec. 

All  the  values  quoted  above  have  been  for  a distance  of  8 cm  downstream  of 
the  beam  and  are  representative  of  the  average  gas  properties  within  the  field  of 
view.  There  can,  however,  be  a significant  variation  in  these  quantities  across 
the  field  of  view.  For  example,  the  values  of  a will  vary  nonlinearly  between 

1.  8 X 10^^  - 4.  8 X 10^^  over  the  range  of  f =4-12  cm.  However,  since  the  | 

fluorescence  measurements  are  integrated  over  the  field  of  view,  these  variations 

will  be  unimportant  as  long  as  beam-produced  species,  such  as  vibrationally  cold 

CO,  do  not  play  a role  in  the  vibrational  relaxation  phenomenology.  It  remains 

then  to  determine  the  local  concentrations  of  beam  produced  species  which  build 

up  during  the  gas  residence  time.  , 

The  gas  flow  within  the  test  chamber  is  somewhat  complicated.  The  test 
chamber  is  cylindrical  with  length  of  16  cm  and  radius  of  10  cm.  The  electron  j 

beam  is  injected  along  the  cylinder  axis;  however,  the  gas  flow  is  not  coaxial.  j 

Gas  is  injected  from  a port  on  the  cylinder  sidewall  situated  5.  5 cm  from  the  beam  i 

origin,  or  alternately,  2.  5 cm  from  the  center  of  the  field  of  view.  There  is  a 
secondary  pumping  source  at  the  beam  origin;  however,  the  primary  pumping 

source  is  near  the  rear  of  the  chamber  on  the  same  cylinder  sidewall  as  the  gas  f 

injection  port.  Thus  gas  flowing  downstream  need  not  remain  on  axis  and,  there-  jj 

fore,  gas  irradiated  by  one  pulse  may  no  longer  be  within  the  radius  of  the  field  ! 

of  view  when  the  beam  is  pulsed  again.  Furthermore,  the  quasistagnant  gas  be- 
tween the  beam  origin  and  the  gas  injection  port  is  not  expected  to  flow  within  the 

field  of  view.  These  statements  cannot  be  firmed  up  without  a detailed  fluid  dynam-  [. 

ical  analysis  which  is  beyond  the  scope  of  the  present  effort.  j 

For  purposes  of  calculation,  it  will  be  assumed  that  the  gas  residence  time  (| 

may  be  defined  in  terms  of  the  time  it  takes  to  traverse  the  field  of  view  at  the 

i! 

average  gas  velocity.  This  latter  quantity  may  be  deduced  to  be  25  cm/sec  from  :i 

i 

the  known  flow  rate  of  5 st.  liters /min,  and  thus  the  time  to  traverse  the  field  of  ; 

view  is  0.  33  sec.  During  this  time  the  beam  would  be  pulsed  13  times  and  there- 
fore, on  the  average,  gas  within  the  field  of  view  will  have  been  irradiated  by  / 

seven  pulses.  The  number  of  CO  molecules  created  in  one  pulse,  over  a length  f, 

is  defined  by  j 

i 
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(60) 


[ 


p(l/p  dE/dx)  I / T 

CO  ^ ^ 

^^tot  e w 

or  for  the  present  case  with  a pulse  time  of  2.4  msec,  current  of  0.  9 mA,  and 
length  of  8,  2 cm 


COjot  = 5 X 10^^  . 

Thus,  the  total  number  of  cold  CO  molecules  within  the  field  of  view  at  one  time 
will  be  3.  5 X 10^^.  Neglecting  diffusion,  one  half  of  these  CO  molecules  would  be 
within  the  gas  volume  defined  by 

12.  2 

V = y*  ’^*‘1/2^  d i ~ 110  cm^  <^2) 

3.  9 


where  Eq.  (52)  has  been  applied.  Thus,  the  average  concentration  of  cold  CO 
within  the  beam  half  current  volume  is  1.  6 X 10^^/cm^.  The  total  beam  irradiated 

volume  defined  by  the  contour  2 r,  is  440  cm^  and  the  average  cold  CO  concen- 

’ 13  3 

tration  in  this  region  would  be  0.  8 X 10  /cm  . 

This  density  of  cold  CO  is  sufficiently  large  so  that  the  reaction 


CO(v)  + CO(0)  -CO(v  - 1)  + CO(l) 


(50) 


could  contribute  a 10  percent  effect  to  the  relaxation  of  individual  vibrational  levels 
of  CO  (see  the  discussion  in  Section  2)  and  is  somewhat  larger  than  desirable. 
However,  given  the  uncertainty  in  the  test  chamber  flow  properties  and  the  neglect 
of  diffusion  effects  (the  characteristic  diffusion  length  being  2 cm  in  0.  3 sec),  the 
estimated  cold  CO  concentration  is  most  probably  an  upper  bound.  A test  of  the 
effect  of  cold  CO  could  be  provided  by  observing  the  relaxation  behavior  of  two 
cases  having  the  same  experimental  conditions,  but  with  beam  currents  differing 
by  a factor  of  ~3.  Such  a test  has  not  been  performed  during  the  present  report- 
ing period. 

The  last  quantity  to  be  determined  is  the  temperature  increase  in  the  gas 
caused  by  beam  heating.  The  average  energy  added  per  cm^  during  the  gas  resi- 
dence time  is  specified  by  the  product  of  the  average  cold  CO  concentration  and 
the  energy  required  to  create  an  ion  pair,  w.  The  average  energy  added  in  the 
beam  irradiated  volume  is  then  3.3  X lO"^  if  err?  which  could  result  in  a 3°K 
increase  in  gas  temperature.  The  temperature  increase  within  the  half  beam 


( 


I 


7\ 
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irradiated  volume,  the  hottest  portion  of  the  gas,  would  be  G^K.  These  minimal 
temperature  increases  should  have  no  effect  on  the  data  analysis. 

3.2  Data  Keduction 

Although  the  predicted  cold  CO  concentration  is  marginal,  the  data  should  be 
sufficient  to  allow  determination  of  the  rate  constants  for  the  process 


CO(v)  + COg  CO(v  - 1)  + COgli/g)  . (3) 

1 

The  raw  .data  is  fluorescent  radiation  arising  from  the  fundamental  vibration/ 
j rotation  band  of  CO.  The  first  step  in  the  data  analysis  is  to  deduce  the  time 

dependent  vibrational  population  histories  from  these  spectra.  Since  the  vibra- 
tional distributions  are  nonBoltzman,  this  step  requires  an  iterative  comparison 
f between  the  data  and  artificially  generated  spectra.  Such  comparisons  can  be 

greatly  simplified  through  the  use  of  nonlinear  regression  analysis,  which  can 
essentially  provide  the  equivalent  of  a least  squares  curve  fit  to  the  data. 

This  technique  involves  the  generation  of  a set  of  "basis  functions" 
where  v is  wavenumber,  which  are  essentially  the  spectra  arising  from  the  individ- 
ual vibrational  transitions,  v — v - 1.  These  spectra,  having  been  properly  con- 
voluted over  the  optical  detection  system  slit  function,  may  be  related  to  the 
observed  intensity  at  any  wavenumber  by  the  relation 


V 

where  the  quantities  n^  are  proportional  to  the  CO  population  in  level  v.  Given  the 
data  Uv)  and  the  functions  defined  over  some  wavenumber  interval,  the 

quantities  n^  may  be  determined  by  requiring  that  the  square  of  the  differences 
(Ip(i/)  - I(i^))  be  minimized.  A computer  code  has  been  developed  to  perform  this 
task  which  takes  the  data  and  basis  functions  as  input  and  provides  the  relative 
vibrational  populations  as  output. 

A computer  code  also  was  developed  to  generate  the  basis  functions,  which 
are  a function  of  both  temperature  and  the  system  slit  function.  The  intensity  of 

an  individual  vibrational/ rotational  transition  of  a diatomic  molecule  with  zero 

37 

angular  momentum,  such  as  CO,  is  defined  by  the  relationship 


37.  Herzberg,  G.  (1951)  Molecular  Spectra  and  Molecular  .Structure  I.  .Spectra 
of  Diatomic  Molecules.  2nd  edition,  D.  Van  Nostrand  Co. , New  York. 
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1 


where  the  vibration/ rotation  interaction  factor  has  been  taken  to  be  unity.  In  Eq, 

((>4)  h is  Planck's  constant,  c is  the  speed  of  light,  ^ j the  H6nl  - London 

factor.  Ay  ^ j is  the  Einstein  coefficient,  J is  the  rotational  quantum  number, 

,,  , , , , is  the  line  center  wavenumber  of  the  transition,  r , is  the 

V — v-l,J  -Jil  v—v-i 

wavenumber  of  the  band  origin,  and  N -is  the  population  of  the  radiating  state. 

37  V,  d 

The  quantity  N - may  be  rewritten 

V,  d 

hcB 

N T = N — (2J  + 1)  exp  (-F(v,  J)  hc/kT)  , (65) 

V,  J V K i 

where  is  the  rotational  constant  for  vibrational  level  v,  defined  as 

B - B -a  (v  + 1/2)  (66) 

V e e 

and 

F(v,  J)  = J (J  + 1)  By  - (J  + 1)^  . (67) 

The  various  line  and  band  origin  wavenumbers  in  Eq.  (64)  may  be  determined  from 

37 

the  definition  of  the  energy  of  a given  v,  J state,  that  is, 

Ey  J = Wg(v  + 1/2)  - Wg  Xg(v  + 1/2)^  + Wg  y^iv  + 1/2)^ 

+ B J(J  + 1)  - D (J  + 1)^  . (68) 

V e 

37 

Finally,  the  HOnl  - London  factors  are  defined  as 

S,  T 1 “ J R branch 

J *J  ” 1 

(69) 

S,  T . 1 J + 1 P branch 

The  intensity  of  a given  v,  J-v  - 1,  J ± I transition  per  vibrationally  excited 
molecule,  that  is,  may  be  defined  from  Eqs.  (64) -(69)  once  the 

Einstein  coefficients  are  defined.  Predictions  for  the  Einstein  coefficients  for  the 
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CO  fundamental  and  first  overtone  bands  as  performed  by  Young  and  Eachus  and 

O Q 

by  Docken  are  shown  in  Figure  16.  As  can  be  seen,  the  difference  between  the 

39 

two  predictions  is  small  and  the  results  of  Docken  have  been  used  in  the  present 
analysis.  It  is  of  interest  to  note  that  the  results  of  Figure  16  are  closely  approx- 
imated by  the  relationship 


^v^v-1 


w 

- 2 V w X \ 

e 

e e \ 

w 

- 2 w X / 

e 

eel 

(70) 


The  line  intensities  were  evaluated  utilizing  the  most  accurate  definition  of  the 

40  -1 

spectroscopic  constants  of  CO  available.  The  values  used,  in  cm  , are 

w =2169.812,  w X 13,289,  w y 1.059  X 10"^,  R 1.  93127, 
e ee  ’ee_,.  'e 

o - 1.  7459  X 10~2  and  D 5.  6268  X lO’  ’.  The  intensities  were  integrated  over 
e e 

the  system  slit  function  g(r,  r')  in  order  to  define  the  basis  functions  f ^(i  )•  -Since 
the  transitions  are  at  discrete  points,  the  integral  becomes  a sum  defined  by 


1 


^'1  =ir  2 e<^'‘'v-v-i. 


. .v  - V - 1 , V 

j-j  ±r  ti  '“'v  - V-  i.j  -j  ±i' 


(71) 


The  slit  function  for  the  system  interferometer  was  a sine  function  defined  by 


0.82Ar 

^ nv  - i/') 


M(r  - I ')  \ 

r9':82Ar  > 


(72) 


where  Ar  is  the  data  spectral  resolution  of  10  cm 

The  basis  functions  ?y(c)  were  generated  for  a temperature  of  300“k  and  used 
to  fit  the  data.  Typical  data  for  the  case  ^002  torr,  1 0.9  niA  are  shown 

in  Figure  17.  These  spectra  were  taken  after  beam  termination  at  times  of  3,  4, 

5 and  6 msec  in  order  of  decreasing  intensity,  respectively,  h'or  purposes  of 
reference,  it  is  noted  that  the  v 1-0  band  origin  is  at  2143  cm  ■*  and  the 
V = 10  - 9 band  origin  is  at  1900  cm  It  can  be  seen  that  the  data  is  quite  noisy. 
The  signal  to  noise  ratio  is  particularly  bad  below  1850  cm  \ corresponding  to 
the  spectral  region  of  decreasing  detector  response.  Thus,  although  spectral  data 
is  available  over  the  band  region  of  the  first  14  vibrational  levels  of  CO,  only 
11-12  vibrational  levels  are  sufficiently  well  defined  to  allow  evaluation  of  their 
population  densities. 

38,  Young,  L. A,,  and  Eachus,  W, J.  (1966)  J.  Chem.  Phys.  44 :4 195. 

39,  Docking,  K.  K.  (1976)  Unpublished  results. 

40,  Schiffner,  G.,  and  Klement,  E,  (1975)  Appl.  Phys.  6:199. 
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Figure  IG.  Predicted  Einstein  Coefficients 
for  the  Fundamental  and  First  Overtone 
Bands  of  CO  (-ref.  38, ref.  39) 


Figure  17.  Normalized  Intensity  at  Times  oft  - 3,  4, 
5 and  G msec  in  Order  of  Decreasing  Intensity,  He- 
spectively.  P,-.,.,  1/4  torr,  I ' 0.  9 m.\ 
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Curve  fitrf  of  fhe  data  were  firjt  performed  in  .-jtep.-j  of  0,  4 msec  between  the 
time  interval  of  2.  6-8.  2 msec.  The  beam  is  off  during  this  time  scale  and,  there- 
fore, the  data  may  be  interpreted  in  terms  of  reaction  <3).  Typical  comparison 
between  the  data  and  curve  fits  at  different  times  are  shown  in  Figures  18-23.  The 
sequence  of  curves  at  the  bottom  of  each  figure  are  the  individual  basis  functions 
and  indicate  the  contribution  of  each  vibrational  level  to  the  total  intensity.  As  can 
be  seen,  the  comparison  between  data  and  prediction  is  quite  reasonable.  In  par- 
ticular, the  predicted  spectra  in  the  region  of  the  lowest  vibrational  levels  is  in 
very  close  conjunction  with  the  data.  This  is  particularly  true  in  Figures  18  and 
19  which  exhibit  the  largest  densities  of  the  v 1 and  2 levels.  The  intensity  slope 
in  the  region  2120-2220  cm  ^ is  due  predominantly  to  the  R branches  of  v 2 -»  1 
and  V 1-0  transitions,  and  the  good  agreement  between  data  and  prediction  con- 
firm that  the  CO  rotational  temperature  is  ®300"k. 

The  comparison  between  data  and  prediction  is  not  as  good  in  the  lower  wave- 
number  region,  £1900  cm  \ although  the  curve  fit  envelopes  the  data  fluctuations 
quite  well.  Severe  noise  fluctuations  are  evident  in  all  the  data  below  = 1840  cm 
Any  errors  in  the  magnitude  of  the  basis  functions  in  this  wavenumber  region  tend 
to  propagate  to  higher  wavenumbers  because  of  the  large  overlap  between  the  basis 
functions.  An  attempt  to  improve  the  fits  by  truncating  the  data  at  1840  cm  ^ 
failed.  The  reason  for  this  is  that  the  magnitude  of  the  basis  functions  in  this 
region  are  overestimated  by  the  fit,  because  there  is  an  actual  contribution  to  the 
intensity  in  this  region  due  to  higher  vibrational  levels.  The  effect  of  this  over- 
estimate propagates  to  higher  wavenumbers  causing  an  underestimate  of  the  magni- 
tude of  the  basis  functions  centered  in  that  region.  It  has  been  found  that  the  curve 
fitting  technique  is  quite  accurate  if  the  band  spectra  is  closed,  that  is,  if  the 
intensity  decays  to  zero  at  both  high  and  low  wavenumber  edges  of  the  band.  Un- 
fortunately, the  noise  level  is  sufficiently  high  in  the  low  wavenumber  range  of 
interest  to  preclude  a reasonable  extrapolation  to  zero  intensity.  In  any  event,  it 
is  felt  that  the  population  densities  of  the  first  12  vibrational  levels  of  CO  can  be 
well  determined  from  the  curve  fits  and  that  the  densities  of  higher  levels  are 
minimal. 

It  can  be  seen  from  examination  of  Figures  18-23  that  the  high  and  low  vibra- 
tional levels  of  CO  decay  more  rapidly  than  those  at  the  center,  that  is.  with 
increasing  time  the  spectra  collapses  towards  its  center.  This  behavior  is  mani- 
fested in  Figure  24  where  the  deduced  relative  vibrational  population  distribution 
is  shown  at  three  different  times.  Note  that  there  is  a total  population  inversion 
in  the  region  of  levels  5-9  and  that  the  populations  of  levels  7-9  vary  only  slightly. 

On  the  other  hand,  significant  decreases  in  the  population  of  higher  and  lower  levels 
can  be  observed.  The  time  dependent  histories  of  the  relative  vibrational  populations 
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Figure  24.  Deduced  Relative  Vibrational  Popula-  / 

tion  Distribution  at  Three  Different  Tiines. 

- 1/4  torr,  I = 0.  9 mA 

are  shown  in  Figures  25  and  2().  These  densities  were  deduced  at  time  intervals 
of  0.4  msec  and  have  been  arbitrarily  connected  by  straight  lines.  As  can  be  seen 
in  Figure  25,  the  first  five  levels  of  CO  relax  smoothly,  exponentially  to  first  order, 
with  the  fluctuations  at  the  longer  times  due  most  probably  to  the  lower  signal  to 
noise  ratio  (S/N)  as  observed  in  Figures  22  and  23.  The  time  histories  for  vibra- 
tional levels  fi-12,  as  shown  in  Figure  20,  are  more  complicated.  Considering 
first  the  highest  vibrational  levels,  it  can  be  seen  that  the  observed  vihratifinal 
relaxation  times  decrease  with  increasing  vibrational  level  for  levels  v 10-12. 

The  results  for  level  12  are  marginal  since  the  band  origin  for  the  v 12-11 
transition  is  at  1840  cm  , or  about  the  point  where  the  signal  to  noise  fluctuations 
become  quite  large,  and  it  is  felt  that  the  populations  for  this  level,  for  limes 
>6  msec,  are  misleading  and  should  be  ignored.  The  behavior  of  level  0 is  quite 
distinct  from  that  of  levels  10-12,  and  in  fact  the  population  of  this  level  actually 
increases  over  the  time  interval  2.  0-5.0  msec.  Furthcrnu'rc,  the  concentration 
of  levels  7 and  8 remain  relatively  constant  over  the  time  interval  shown  (and 
beyond)  although  the  population  of  level  0 drops  by  a factor  of  two. 


Figure  25.  Deduced  Time  Histories  of  the  Relative 
Populations  of  the  First  Five  l.evels  of  CO. 
torr,  1 = 0.9  mA  2 
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Figure  2(1.  Deduced  Time  Histories  of  the  Relative 
Populations  of  l.evels  v (>-12  of  CXI.  Pfxis  ^ ^ 
tore,  I 0.  9 mA 


Time,  msec 


These  observations  imply  that  levels  7-9  act  as  a bottleneck  for  efficient 
relaxation  of  vibrationally  excited  CO  in  this  system.  This  behavior  can  best  be 
interpreted  through  an  analysis  of  the  kinetic  mechanisms  controlling  the  relaxa- 
tion phenomena.  In  the  absence  of  intramode  vibrational  exchange  collisions  - 
that  is,  processes  (4)  - the  master  equations  for  relaxation  of  vibrationally  excited 
CO  may  be  written  as 


d = ^v+1  %+l  • ^vV  ^ ^v+2-v  "v+2  - Av-v-2  "v 

for  V = 1 - V . In  Eq.  (73),  n is  the  relative  population  of  level  v,  A „ 

max  V v-r^-*v 

is  the  Einstein  coefficient  for  first  overtone  radiative  transitions,  as  shown  in 
Figure  16,  and  is  defined  as 

Rv  = VCO2)  + . (74) 

where  k is  the  rate  constant  for  reaction  (3),  (CO,)  is  the  CO,  number  density  in 
Vo  i 

particles/cm  , and  v 1 Einstein  coefficient  for  fundamental  band  radia- 

tive transitions.  The  unknown  quantities  in  Eq.  (73)  are  the  k^'s  and  the  derivatives, 
dn  /dt.  In  theory,  the  derivatives  could  be  determined  from  the  predicted  time 
variation  in  n^;  however,  in  practice  this  technique  is  not  tractable  because  of  the 
uncertainties  in  the  predicted  n^'s.  Therefore,  Eq.  (73)  has  been  used  in  its  time 
integral  form  which  may  be  written 


nvitg)  - "^(*1)  - ^v+1  Rv+1 


a R + A , , 
V V v+2- 


-2  ®v+2  ■ ^v- 


v-2  ^v 


(75) 


where 


n^dt 


(76) 


Since  the  values  of  a^  can  be  determined  directly  from  the  data,  the  quantities  R^ 
may  be  deduced  from  Eq.  (75). 

This  determination  can  be  simplified  by  the  fact  that 


('^v*‘2'  ■ "v^4^'  ■ '^1^1  ■ ■^2- 


0 ®2 


v=l 


(77) 
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which  can  be  shown  directly  from  Eq.  (75).  Equation  (77)  simply  states  that  the 
decrease  in  the  population  of  vibrationally  excited  CO  can  only  occur  through  vibra- 
tional deactivation  of  level  1 and  first  overtone  radiation  from  level  2 (radiation 
from  higher  overtones  being  negligible).  Once  the  quantity  Rj  is  defined,  the 
remaining  R^'s  can  be  determined  through  the  relationship 

n-1 

2 = -^1«1  -^2,0^2  ^^n^n^  Vl^n-1 

v=l 

The  complication  in  using  Eqs.  (77)  and  (78)  is  that  the  system  of  equations  is  not 
closed;  that  is,  relaxation  data  is  only  available  for  levels  v < 12  and  level  12 
itself  is  not  defined  very  accurately.  Although  the  behavior  of  levels  higher  than 
12  will  not  effect  the  determination  of  Rj^,  since 

it  will  affect  the  determination  of  R^  for  large  v,  inasmuch  as  these  will  be  a small 
difference  between  two  large  numbers.  Therefore,  the  values  for  R^  determined 
at  higher  levels,  v > 9,  must  be  considered  as  lower  bounds. 

The  values  of  R as  determined  from  the  data  of  Figures  25  and  26  and 
Eqs.  (77)  and  (78),  are  shown  in  Figure  27.  The  values  shown  for  Rj2  has  been 
deduced  directly  from  the  observed  decay  of  that  level.  Also  shown  for  compari- 
son are  the  values  of  the  Einstein  coefficient  A The  R 's  shown  have  been 

evaluated  for  three  different  time  intervals,  as  shown.  The  largest  time  interval 
is  to  be  preferred  inasmuch  as  this  case  will  minimize  the  effect  of  uncertainties 
in  the  data.  Similar  results  are  obtained  if  tj  is  taken  to  be  3 msec.  It  would 
appear  from  the  data  of  Figure  27  that  one  unique  set  of  R^'s  is  consistent  with 
the  relaxation  observations,  (The  differences  in  the  R^'s  determined  at  different 
time  intervals  are  only  pronounced  at  large  v and  these  are  most  affected  by  small 
variations,  as  discussed  earlier, ) 

It  is  also  evident  from  Figure  27  that  the  rate  constants  k^  cannot  be  deter- 
mined with  any  accuracy  for  levels  v = 8,  9.  This  is  because  R^  * A^  ^ j for 
these  levels  and  the  rate  constant  k^  is  proportional  to  the  difference  between  R^ 
and  A^  (see  Eq.  (74)).  Furthermore,  the  R^'s  for  levels  v = 7 and  10-12  are 
all  within  a factor  of  two  of  the  respective  values  of  A^  Coupling  this  obser- 

vation with  the  fact  that  the  values  shown  for  R^  at  the  higher  levels  are  lower 
bounds  — because  the  effect  of  levels  higher  than  v ^ 12  has  not  been  included  in  the 
analysis  — it  is  clear,  that  the  uncertainty  in  the  rate  constants  for  these  levels  will 
also  be  high. 
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Figure  27.  CC)(v)  Vibrational  Helaxatiun  KatH  i- 
P(-'02  ^ torr,  I 0.  !•  niA.  .A^.  ' 'rr 

bound  from  data 


The  nominal  pressure  for  this  case  is  1 4 torr;  however,  it  is  uncertain 

to  « 50  percent.  The  exact  CO2  pressure  has  been  determined  from  the  data  by 

setting  kj  equal  to  its  known  value  as  shown  in  I igure  10.  The  pressure 

determined  in  this  manner  was  1/3  torr.  This  pressure  and  the  measured  values 

of  were  then  used  to  evaluate  k^  through  use  of  F.q.  (74).  These  results  are 

shown  in  Figure  28  for  the  three  time  intervals  used.  Also  shown  are  the  results 
24 

of  Hancock  and  Smith.  As  can  be  seen,  there  is  a large  scatter  in  the  rate  con- 
stants V a 7,  while  those  for  v < 7 appear  well  defined.  The  rate  constants  for 
v < 7 are  considered  to  be  reasonably  accurate;  however,  they  tend  to  be  a factor 
of  2-3  larger  than  the  rate  constants  as  determined  by  Hancock  and  Smith.  On 
the  other  hand,  the  rate  constants  for  v ■ 7 are  smaller  than  those  of  Hancock  and 
Smith,  as  would  be  expected  inasmuch  as  they  are  lower  bounds.  The  cause  of 
the  difference  between  the  rate  constants  for  v 4-()  as  determined  in  this  analysis 
and  Hancock  and  Smith  is  not  understood.  However,  the  present  analysis  is  most 
accurate  at  low  vibrational  levels  whereas  that  of  Hancock  and  Smith  is  most  accu- 
rate at  high  vibrational  levels,  f'urthermore,  although  if  is  recognised  (hat  file 
presence  of  cold  CO  could  cause  an  overestimate  of  the  rate  constants,  it  is  unlikely 


I I CO2  - CC)(V  -1)4  C02(i;^)  o.  □.  0 ••  present 
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that  this  effect  could  cause  differences  of  factors  of  2-3.  Nonetheless,  it  is  clear 
that  additional  cases  must  be  reduced  in  order  to  check  the  consistency  of  the  pres- 
ent results. 

Two  internal  consistency  checks  have  been  performed.  The  first  was  a com- 
parison of  the  measured  vibrational  population  time  histories  with  those  predicted 
by  solution  of  the  differential  Eq.  (73),  using  the  deduced  values  of  R^.  Tlie  solu- 
tion of  these  equations  is  also  complicated  by  closure  effects  inasmuch  as  the 
population  of  levels  v 13  and  14  — which  arc  unknown  — are  needed  for  the  evalua- 
tion of  the  differential  equation  for  the  population  of  v 12.  This  complication  lias 
been  circumvented  by  solving  the  equations  for  v 1-11  taking  the  population  his- 
tory of  V 12  as  known  from  the  data  (the  effect  of  v 13  first  overtone  decay  to 
V 11  is  neglected).  I'urthermore,  the  population  of  v 12  was  taken  to  continue 
to  decay  exponentially  beyond  t fi.  fi  msec  inasmuch  as  the  cui've  fit  results  are 
suspect  at  the  longer  times.  (The  value  of  R,.,,  as  determined  by  this  exponential 
decay,  is  430  sec  .) 

The  comparison  between  predicted  and  measured  time  histories  for  the  first 
five  vibrational  levels  is  shown  in  Figure  29  and  can  be  seen  to  be  <)uite  good.  The 


the  assumed  time  history  of  level  v = 12  used  in  the  calculations.  Once  again  the 
comparison  between  prediction  and  data  is  excellent  for  levels  v = 6-9,  with  all 
the  salient  features  of  the  data  being  exhibited  in  the  predictions.  The  predicted 
time  histories  for  levels  v = 10-11  decay  somewhat  more  slowly  than  the  data 
implying  that  and  Rj2  are  somewhat  higher  than  deduced  in  the  prior  analysis. 
The  overall  good  comparison  with  the  time  depem'.ent  observations  suggests  that 
the  kinetic  mechanism  used  in  interpreting  the  data  was  appropriate. 

The  second  check  involved  the  total  vibrational  energy  decay  in  CO.  The  total 
radiation  from  the  fundamental  band  is,  to  first  order,  linearly  proportional  to  the 
vibrational  energy  content  of  the  molecule.  Therefore  the  experimental  spectra 
were  integrated  over  the  spectral  range  of  1760-2240  cm  \ and  this  time  dependent 
fundamental  band  radiation  was  compared  to  the  predicted  decay  in  vibrational 
energy  as  determined  from  the  computer  solution  of  Eq.  (73).  This  comparison  is 
shown  in  Figure  31  and  is  seen  to  be  quite  reasonable.  The  somewhat  more  rapid 
I decay  observed  in  the  data  at  long  times  is  expected  inasmuch  as  the  decay  rates 

for  higher  levels  were  underestimated.  This  good  comparison  suggests  that  no 
serious  errors  developed  when  the  experimental  spectra  were  converted  to  vibra- 
tional populations. 

Although  it  is  not  clear  why  the  present  measurements  differ  from  those  of 
24 

Hancock  and  Smith,  the  total  single  quantum  deactivation  rates  occurring  in  the 
existing  gas  mix  have  been  measured.  These  may  then  be  used  to  analyze  the 
fluorescence  data  taken  while  the  beam  was  on,  t = 0-2.  4 msec,  in  order  to  deter- 
mine the  production  rates  for  vibrationally  excited  CO.  Typical  beam-on  spectra 
along  with  curve  fit  artificial  spectra  are  shown  in  Figures  32-34.  These  data  and 
fits  are  very  similar  to  those  discussed  earlier.  The  spectra  at  the  lowest  times 
are  quite  noisy  and  the  curve  fit  predictions  for  t S 0.  6 msec  must  be  considered 
as  crude  estimates.  These  curve  fits  were  performed  in  steps  of  0.  2 msec  for 
times  between  0.4-2.  4 msec,  and  the  resulting  relative  vibrational  populations 
are  shown  in  Figures  35  and  36.  Also  shown  for  comparison  are  the  previously 
determined  vibrational  populations  shortly  after  beam  termination.  As  can  be 
seen  in  Figure  35,  the  populations  of  levels  1-5  increase  monotonically  while  the 
beam  is  on,  with  populations  decreasing  with  increasing  vibrational  level.  The 
transition  to  vibrational  decay,  upon  beam  termination,  occurs  smoothly.  I’erhaps 
the  most  striking  behavior  is  seen  in  the  population  histories  of  levels  6-12,  shown 
in  Figure  36.  Here  it  can  be  seen  that  the  population  for  level  v ' 9 is  less  than 
that  for  V - 10-12  at  early  times.  However,  near  beam  termination,  the  population 
of  V - 9 continues  to  increase  while  those  of  the  higher  levels  go  through  a 
maximum. 
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Figure  31.  Comparison  of  Measured  and  Predicted  CO  Vibrational  Knergy  Decay 


These  population  histories  may  be  analyzed  in  a manner  similar  to  that  used 
earlier  for  the  vibrational  decay  data.  The  governing  set  of  differential  equations 
are 


dn  /dt  R (v)  t R n 

\t’  V4- 1 V 


t A , n . „ - A 
v‘2-v  v+2  V- 


^ n 

v-2  V 


where  R (vl  is  the  relative  rate  for  creation  of  CO  molecules  in  vibrational  level  v 
e 

(this  rate  is  only  relative  inasmuch  as  the  n^  are  relative  vibrational  populations). 
Once  again  the  integral  form  of  Eq.  (HO)  is  used,  that  is, 


n (t>  R (v)  t * a R - a R 
V e v+1  v^l  V V 
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Figure  32.  Synthetic  and  Measured  Spectra  at  t 0.  (>  msec 
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Figure  33.  Syntlietic  and  Measured  Spectra  at  t 1 . (i  msec 
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Figure  34.  Synthetic  and  Measured  Spectra  at  t = 2.  4 msec 

The  results  for  the  time  interval  of  0-2.4  msec  are  shown  in  Figure  37.  As  can 
be  seen  vibrationally  excited  CO  is  created  predominantly  in  the  lower  vibrational 
levels;  however,  there  is  a secondary  peak  in  the  rate  distribution  at  level  v = 6 
and  the  creation  rate  for  v = 9 is  below  those  for  neighboring  levels.  The  rate  for 
creation  of  v * 12  cannot  be  accurately  deduced  from  the  data,  given  the  lack  of 
information  on  the  populations  of  v - 13,  14  and  the  uncertainty  in  Rj2- 

As  a check  on  the  analysis,  the  rates  R^iv)  can  be  used  in  conjunction  with 
Eq.  (80)  to  compute  the  relative  vibrational  population  time  histories  while  the 
beam  is  on.  Such  calculations  are  shown,  in  comparison  to  the  data,  in  Figures  38 
and  39.  The  computed  values  for  the  first  five  vibrational  levels  are  shown  in 
Figure  38  and  can  be  seen  to  be  in  excellent  agreement  with  the  data  at  times  >0.  8 
msec.  For  times  below  0.8  msec,  the  data  for  levels  v - 3-5  falls  off  more 
sharply  than  calculation;  however,  it  is  to  be  remembered  that  the  data  is  highly 
inaccurate  at  short  times  (for  example  see  F'igure  32).  The  comparison  for  levels 
V 8-11  is  shown  in  Figure  39  and  once  again  the  agreement  between  data  and  cal- 
culation is  quite  reasonable.  (The  Eqs.  (80)  were  solved  by  taking  the  population 
of  level  11  fixed  by  the  data,  with  extrapolation  to  zero  time  as  shown  in  Figure  39. ) 
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Figure  37.  Relative  Rates  for  Creation  of  Vibra- 
tionally  Excited  CO.  ^002  " torr,  1 = 0.9  mA 
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Figure  38.  Comparison  of  Measured  and  Computed 
Vibrational  Population  Time  Histories  for  v = l-.'j. 
— data, prediction.  Beamon 
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Figure  39.  Comparison  of  Measured  and  Computed 
Vibrational  Population  Time  Histories  for  v = 6-11. 
— data, prediction.  Beam  on 


3.3  Data  Interpretation 

One  data  set  for  CO  has  been  completely  reduced  and  a set  of  rates  for  crea- 
tion and  deactivation  of  vibrationally  excited  CO  have  been  determined.  It  would 
be  valuable  to  reduce  additional  data  sets  in  order  to  ascertain  whether  the  present 
results  would  be  consistent  with  observations  taken  under  different  experimental 
conditions.  Such  attempts  have  been  unsuccessful  during  the  present  reporting 
period.  In  particular,  strenuous  efforts  were  made  to  develop  CO  vibrational 
population  distributions  from  spectra  taken  at  the  conditions  of  I = 3.  0 mA  and 
=1/4  and  1 torr.  These  efforts  met  with  failure  primarily  because  of  the 
low  signal  to  noise  ratio  of  the  spectra.  Nonetheless,  several  salient  features 
could  be  observed  in  the  data  base;  (a)  at  both  PCO2  " ^ torr  the  total 

fundamental  band  radiation  scaled  linearly  with  beam  current  for  current  of  0.  9 
and  0.  3 mA,  (b)  the  deactivation  rates  for  high  and  low  vibrational  levels  increased 
with  increasing  CO2  concentration  between  l^C02 ' 0-  torr,  and  (c)  the  peak 
CO  intensity  at  beam  termination,  occurring  at  * v = 7-9,  remained  relatively  con- 
stant over  the  CO2  pressure  range  of  0.  1-1  torr.  These  observations  are  all  con- 
sistent with  the  deactivation /excitation  rates  and  mechanisms  deduced  in  the  pres- 
ent analysis. 


J 


1 
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It  would  be  of  value  to  specify  the  absolute  CO  concentrations  and  excitation 
rates  in  order  to  determine  what  portion  of  the  electron  energy  is  ultimately  con 
verted  into  CO  vibrational  energy.  The  absolute  CO  vibrational  densities  within 
the  field  of  view,  may  be  related  to  the  total  band  intensity,  I.p,  specified  in 
Watts/unit  area/sterradian  by  the  relationship 


SNA  , hi/  , 
V V-»V“1  V- 


v-1 


4 irA, 


(82) 


where  A^  is  the  area  of  the  field  of  view.  The  are  in  turn  linearly  proportional 
to  the  relative  CO  number  densities  discussed  above,  that  is. 


N^  = An^ 


(83) 


and  thus,  given  the  relative  vibrational  populations  and  the  total  CO  band  intensity, 
the  proportionality  constant  A may  be  readily  determined.  The  total  creation  rate, 
within  the  field  of  view,  for  vtbrationally  f xcited  CO  may  then  be  specified  by 


CO*  = A S Rg(v)  . (84) 

Unfortunately,  due  to  complications  involving  nonuniform  detector  response 
across  the  field  of  view,  an  absolute  calibration  cannot  yet  be  provided.  However, 
based  upon  the  available  information,  it  can  be  estimated  to  within  an  order  of 
magnitude  that  there  is  one  vibratlonally  excited  CO  molecule  created  per  ion  pair, 
that  is. 


CO*  ..  = 0.9pVjdE/dX)l.  ^ 9 ^ ,917  (85) 

where  i is  the  diameter  of  the  field  of  view  and  0.  9 is  the  fraction  of  beam  elec- 
trons which  fall  within  the  field  of  view  over  this  length.  Furthermore,  the  average 
number  of  quanta,  v,  created  per  excited  CO  molecule  may  be  defined  by 


S vR^(v) 
" S Rg(v) 


3.  5 


(86) 


Equations  (85)  and  (86)  may  be  used  to  estimate  that  approximately  3.  5 percent  of 
the  beam  energy  deposited  in  the  gas  is  converted  to  CO  vibrational  energy.  It  is 
to  be  emphasized  that  this  estimate  is  only  accurate  to  within  an  order  of  magnitude. 
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The  actual  kinetic  mechanism  by  which  this  energy  transfer  occurs  is  not  well 
defined.  Reaction  (26)  is  exothermic  by  8.32  eV,  assuming  is  in  its  ground 

vibrational  state.  This  is  sufficiently  energetic  so  that  the  following  reaction 
paths  could  occur 

C02^  + e - CO(A  V ^ 2)  + O (87a) 

COg’*’  + e - CO  (a  ^TT,  V > 0)  + O AH  > -2.3  eV  (87b) 

COg^  + e - CO  (x  V > 0)  + O AH  > 8. 32  eV  (87c) 

Gutcheck  and  Zipf^^  have  studied  the  recombination  of  CO,^  ions  produced  via 
charge  exchange  from  Ar  , as  in  the  present  system,  and  found  that  ® 5 percent 
of  the  recombining  molecules  formed  the  CO  (A  V)  state.  It  was  also  found  that 
the  A state  was  formed  in  levels  v = 0-3  in  approximately  equal  concentrations. 
(Note,  V = 3 was  formed  because  the  recombining  ions  were  vibrationally  excited. 

It  is  anticipated  that  this  would  not  be  the  case  in  the  present  system. ) The  A 

so  created  would  radiatively  decay  to  vibrationally  excited  ground  state  CO  with  a 

“8  4 1 

lifetime  of  ~10  sec.  The  relative  CO(v)  creation  rates  resulting  from 

CO(A  ^tt)  -•  CO(X  ^S)  transitions  from  equal  populations  of  CO(A  ^ir,  v = 0-2)  and 

CO(A  V = 0-3),  have  been  estimated  from  Franck-Condon  factors  presented 

40 

in  Schiffner  and  Klement  and  are  shown  in  Figure  40  in  comparison  with  t^e 
measured  excitation  rates.  As  can  be  seen  CO(v)  creation  by  this  process  would 
exhibit  a less  steep  dependence  with  vibrational  level  than  observed.  Of  course, 
the  existence  of  equal  populations  in  the  CO  (A  ^tt)  vibrational  levels  can  only  be 

conjectured  at  present. 

3 + 

CO(a  it)  formation  has  also  been  observed  in  CO,  recombination.  Wauchop 

42  + ^ 3 

and  Broida  created  CO,  through  Penning  ionization  collisions  with  He(2  S) 

“ + 

metastables  and  found  that  55  percent  of  the  CO,  formed  in  this  manner  recom- 
3 ^ 

bined  to  form  CO(a  n),  (Again,  extrapolation  of  these  results  to  the  present  con- 
ditions is  questionable  inasmuch  as  the  ions  produced  in  the  experiment  of  Wauchop 
42 

and  Broida  could  have  been  both  electronically  and  vibrationally  excited. ) There 

3 

is  sufficient  exothermicity  in  reaction  (87b)  to  produce  CO(a  ir)  in  vibrational 

states  as  high  as  v = 11;  however,  states  higher  than  v = 4 will  rapidly  convert  to 
3 3 

the  triplet  states  a 2,  d A,  etc.,  through  radiationless  curve  crossings,  and 

3 

these  states  will  then  radiatively  cascade  back  to  the  a it  state.  The  radiative 


41.  Krupenie,  P.  H.  (1966)  The  Band  Spectrum  of  Carbon  Monoxide,  NSRDS-NBS5. 

42,  Wauchop,  T.S.,  and  Broida,  H.  P.  (1972)  J.  Chem.  Phys.  ^:330. 
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Figure  40,  Estimated  Relative  CO(v)  Creation  Rates 
Due  to  CO(A  ^ir)  -*  CO(X  ^2)  Transitions.  6,  meas- 
ured creation  rates  Rp(v),  — CO(A  ^jr,  v = 0-2), 

---  CO(A  Itt,  V = 0-3) 


3 

lifetime  for  the  a t state  is  several  msec  (Wauchop  and  Broida)  but  it  can  be 
rapidly  quenched  by  collisions  with  COg,  that  is, 


CO(a  V < 4)  + CO,  -*  2CO(X  ^2  ) + O 


AH  -:2  - 1.  1 eV 

V 


(88) 


-1  I '(  4?  41 

Process  (88)  has  a rate  constant  of  1.7-3. 8 X 10  cm  /sec  ’ and  thus  for 
Pco  =1/4  torr,  will  have  a characteristic  quenching  time  of  ~ 5 osec.  The 
exothermicity  of  reaction  (88)  is  sufficient  to  produce  ground  state  CO  in  levels 
V ^ 4,  and  thus  this  reaction  alone  is  not  sufficient  to  explain  the  present  observa- 
tions. It  is  Interesting  to  note,  however,  that  the  measured  CO(v)  creation  rates 
plateau  at  v = 4-5,  with  the  dominant  production  occurring  in  the  first  four  levels. 
Thus,  it  is  possible  that  both  the  CO(A  ^;r)  and  CO(a  ^t)  states  contribute  to  CO(v) 
creation  in  the  present  experiment.  This  conjecture  could  be  tested  by  monitoring 
the  VUV  radiation  occurring  in  the  gas  mixture  while  the  beam  is  on. 


43.  Slanger,  T.G.,  and  Black,  G.  (1971)  J.  Chem.  Phys.  ^:21G4. 
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In  conclusion,  experimental  studies  (Krupenie,  and  Wauchop  and  Broida  ) 

have  shown  that  = 60  percent  of  recombinations  lead  to  the  formation  of 

excited  electronic  states  of  CO.  Although  the  ions  in  these  studies  were  not  neces- 
sarily in  the  same  state  as  those  in  the  present  work,  it  has  been  shown  that  these 
electronic  states  can  lead  to  the  production  of  vibrationally  excited  ground  state 
CO  either  by  radiative  decay  or  by  collisional  quenching. 

t.  CO.,  n\T\  WAl.VSIS 

Little  has  been  said  so  far  about  the  COg  fluorescence  data.  It  wi'i  be  shown 
that  the  002(1^3)  band  radiation  is  severely  optically  thick  over  the  fvil  range  of 
experimental  conditions.  This  effect  precludes  the  determination  ol  the  absolute 
CO2  concentrations  from  data  taken  with  the  integrating  sphere  in,  inasmuch  as  the 
transmission  path  length  within  this  device  is  not  well  defined.  The  data  discussed 
in  Section  3,  taken  without  the  integrating  sphere  in,  are  amenable  to  reduction  to 
absolute  magnitude  once  the  transmission  losses  are  defined.  In  any  event,  the 
raw  fluorescence  data  may  be  used  to  determine  the  deactivation  mechanisms 
occurring  within  the  gas  after  beam  termination. 

C02(r3)  band  fluorescence  decay  data  have  been  taken  at  nominal  CO2  pressures 
of  0,  0.  1,  0.  25,  0.  5,  and  1.  0 torr.  The  zero  pressure  case  corresponds  to  pure 
Ar  with  a minimal  CO2  concentration  resulting  from  wall  out-gassing,  and  the 
C02(i'3)  radiation  in  this  case  is  optically  thin.  At  all  higher  pressures,  the 
radiation  will  be  optically  thick.  The  radiation  decay  in  all  cases  is  exponential, 
to  first  order,  with  characteristic  decay  time  varying  slightly  with  CO2  pressure. 
The  dominant  deactivation  mechanisms  in  the  gas  are  expected  to  be 

C02(i'3)+Ar  002(3^2)  + Ar  kg^  = 1.  7 X 10'^^  cm^/sec  (89) 

^qn  -14  ? 

002(1/3)  + CO2  - C02(3i/2)  + CO2  kgp  = 1.  1 X 10  cm'’ /sec  (90) 

and 

002(^3)  " CO2  + hr  = 2,3  msec  (91) 

44 

where  the  room  temperature  rate  constants  listed  are  from  Yardley  and  Moore. 

The  characteristic  decay  time  for  may  then  be  defined  as 


44.  Yardley,  J.T.,  and  Moore,  C.  B.  (1967)  J.  Chem.  Phys.  46:4491. 
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(92a) 


' <*^89  *^90  "^R 

-.3  = <‘^89[^^1  ^‘^goicoal)"'  . 

where  the  bracketed  quantities  represent  the  number  densities  of  Ar  and  COg.  The 
comparison  between  predicted  and  measured  003(1^3)  decay  times  is  shown  in 
Figure  41  and  can  be  seen  to  be  quite  good  at  the  lower  CO2  pressures.  Note  that 
the  difference  in  observed  relaxation  times  at  CO2  pressures  of  0.  0 and  0.  1 torr 
illustrate  the  effect  of  optical  thickness.  It  can  be  seen  that  the  measured  relaxa- 
tion time  lags  behind  the  predicted  with  increasing  CO2  pressure.  It  is  anticipated 
that  this  effect  is  due  to  vibrational  pumping  of  €02(1^3)  by  CO,  that  is, 

CO(v)  + CO2  - CO(v  - 1)  + 002(1/3)  (3) 

which  becomes  more  pronounced  with  increasing  CO2  pressure.  This  effect  has 
not  been  modeled  since  the  CO  vibrational  populations  at  the  higher  CO2  pressure 
have  not  been  determined. 

In  general  then,  at  the  lower  CO2  pressures  the  002(1/3)  obser- 

vations provide  validation  of  the  measured  rate  constant  for  reaction  (89).  Since 
the  002(1/3)  deactivation  phenomena  is  well  understood,  the  excitation 

rate  may  be  determined  from  the  data.  It  is  expected  that  collisions  with  second- 
ary electrons  provide  the  relevant  excitation  mechanism,  that  is, 

k 

e + CO2  - e + 002(1/3)  . (93) 

From  Eqs.  (92)  and  (93),  it  can  be  shown  that  the  total  002(1/3)  concentration  with- 
in the  field  of  view  while  the  beam  is  on  may  be  defined  by 

N = R T (1  - exp  (-t/r  ))  , (94) 

I'a  *^3  ^'3 


-1 


CO„ 


CO. 


> 0 


(92b) 


where  R is  the  elobal  excitation  rate  for  reaction  (93).  Thus,  if  N,,  can  be 
i-'S  ' 3 

determined,  the  absolute  CO2  excitation  rate  and  beam  excitation  rate  and  beam 

excitation  efficiency  may  be  specified. 

It  is,  unfortunately,  quite  difficult  to  accurately  define  from  the  experi- 
mental spectra.  The  002(13)  radiation  with  the  chamber  is  dominantly  doppler 
broadened  and  severely  self  absorbed.  After  the  radiation  leaves  the  chamber,  it 
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Figure  41.  Comparison  of  Predicted  and  Observed  COqIi'o)  Relaxation  Times. 
= 9 torr,  ^ - data,  Q " prediction  from  Eq.  (92) 


travels  a length  of  30-45  cm  in  ambient  air  before  striking  the  detector  and  there- 
fore undergoes  additional  absorption  by  Lorentz -broadened  ambient  CO2.  This 
latter  effect  has  been  rather  cavalierly  accounted  for  by  application  of  a low  resol- 
ution calibration  correction.  The  justification  for  this  step  will  be  provided  further 
in  the  text.  The  self  absorption  occurring  within  the  chamber  has  been  modeled  in 
some  detail  as  discussed  below. 

It  was  shown  in  Figure  15  that  the  radiating  gas  will  be  localized  within  the 
field  of  view.  For  purposes  of  calculating  the  transmittance,  it  has  been  assumed 
that  the  vibrationally  hot  gas  is  present  only  along  the  axis  of  the  test  chamber  and 
propagates  a distance  of  11  cm  through  a room  temperature,  vibrationally  equi- 
librated gas  before  leaving  the  chamber. 

The  fluorescent  lines  emitted  by  the  excited  CO,  molecules  will  have  a doppler 

_7  ^ 

halfwidth  at  half  maximum  (IfWHM)  of  9.4  X 10  r (where  u is  the  frequency  of  the 

-3  -1 

transition  in  wavenumbers),  or  bp  = 2.  2 X 10  ' cm  HWHM  for  the  1/3  transition 
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of  COg.  Collisional  broadening  of  the  spectral  line  by  Ar  results  in  a broadening 
of  5.  8 X 10'^  cm'*^  HWHM  for  P.  „ = 9 torr.  CO,  self-broadening  will  add 

C A 1 ^ ^ “ 

1.3  X 10"^  - 1.3  X lO'  cm"  HWHM  to  the  width  of  each  line  depending  on  amount 
present  in  the  experiment  (0.  1-1  torr).  Therefore,  doppler  broadening  is  the 
dominant  mechanism  for  spectral  line  broadening,  and  will  be  considered  first  in 
the  modeling  to  follow.  The  unexcited  part  of  the  gas  mixture  is  taken  to  have  the 
same  temperature  and  pressure  as  the  radiating  molecules,  and  thus  the  absorption 
linewidths  will  be  identical  to  the  emitted  ones. 

A doppler  lineshape  may  be  described  by  the  formula 


where  k is  the  absorption  coefficient  at  frequency  r.  k , the  absorption  coef- 
ficient  at  line  center,  is  defined  by 


k Jin2/7r  (96) 

-Cl  S 


where  S is  the  integrated  linestrength  for  the  single  rotational  transition  in  cm/ 
molecule,  that  is,  ^ ~ J which  contains  the  lineshape  behavior  with  fre- 

quency variation,  is  defined  by 


(97) 


It  is  seen  from  these  equations  that  k = 4 k when  v - v - b_.  For  an 

V 2 Vj.  ° 

Lj 

optically  thin  source,  the  radiation/unit  length  of  an  emitter  with  concentration 
Np  molecules/cc  is 

I = k N„  W/cm^-sr-cm  ^ (98) 

v 1/  W DU 

where  Igg  is  the  blackbody  function.  Attenuation  of  this  radiation  by  doppler 
broadened  lines  over  a distance  L is  given  by 


I ' = k N„I„_  exp  (-KNL) 

1/  V DD 


(99) 
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where  K is  the  total  absorption  coefficient  due  to  all  transitions  (of  CO2  or  other 
species  present)  which  overlap  the  emitted  line,  that  is. 


(100) 


The  fractional  transmission  of  a single  doppler  rotational  line  is  then  given  by 


T 


/ ‘v’ 


dv 


/ 

-flO 


exp  ( -K  N L)  dv 


-00 


(101) 


Determination  of  the  transmission  function  requires  knowledge  of  all  at  fre- 
quency intervals  small  compared  to  the  linewidth  bj^  and  evaluation  of  Eq.  (101)  at 
a large  number  of  frequencies  for  each  rotational  transition.  This  approach  is 
both  time  consuming  and  costly.  If  the  assumption  is  made  that  the  emission  from 
a particular  line  is  only  significantly  absorbed  by  the  same  transition  in  CO2  ihat 
gave  rise  to  the  radiation,  a considerable  simplification  in  Eqs.  (100)  and  (101) 
occurs,  that  is,  KN  = k N^j.  Since  the  blackbody  function  remains  constant  over 
the  molecular  linewidth,  Eq.  (101)  then  reduces  to 


dv 


T 


00 


2 


(102) 


Equation  (102)  may  be  represented  as  a series  expansion  for  small  values  of 


k N„  L (Mitchell  and  Zemansky^^): 


E 


k Np  L 


n=0 


11 ! ^ 


(103) 


+ 1 


For  large  values  of  k L,  the  above  series  involves  small  differences 

between  large  terms,  and  numerical  evaluation  of  the  integral  is  required.  A plot 
of  the  fractional  transmission  given  by  Eq.  (101)  as  a function  k L,  is  pre- 

sented  in  Figure  42,  Note  that  t approaches  0 and  1 for  very  strong  and  very  weak 
lines  respectively.  The  total  fractional  transmission  of  the  entire  spectrum  is  just 
the  sum  of  the  transmission  of  each  rotational  line  in  the  band. 

In  order  to  estimate  an  upper  limit  for  the  effect  of  collisional  broadening  on 
the  spectral  lines,  the  fractional  transmission  of  each  rotational  line  was  calculated 
assuming  the  spectral  linewidth  to  be  due  entirely  to  collision  broadening, 
bj.  = 0.0022  + 0.0007  = 0.0029  cm"^  HWHM. 

The  absorption  coefficient  over  a single  collision  broadened  rotational  line  is 
given  by 


cm  = k 


o (v  - v^)  + b^. 


(104) 


where  k is  the  absorption  coefficient  at  the  center  of  the  homogeneous  line  defined 
. *'0 

by 


(105) 


JTb- 


For  the  case  of  an  optically  thin  Lorentzian  source  attenuated  by  Ixjrentzian 
absorbers,  Eq.  (101)  becomes 


45.  Mitchell,  A.C.G.,  and  Zemansky,  M.W.  (1934)  Kesonance  Radiation  and 
Excited  Atoms.  Cambridge,  London,  p.  101. 
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Fractional  Transmission, 


Figure  42.  Fractional  Transmission  as  a Function  of  Absorbing  Gas  Line  Strength 
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(107) 


With  appropriate  substitutions  in  Eq.  (107),  the  transmission  is  foun<I  to  be  a 
function  of  the  zeroth  order  imaginary  Bessel  function,  I^, 

T = e'P  yp)  (108) 

where  p = may  be  approximated  as  a polynomial  permitting  analy- 

tical solution  for  r at  every  value  of  SNjj  L. 

Doppler  and  collision  broadening  are  entirely  independent  processes.  To 
describe  a lineshape  in  which  both  significantly  contribute,  each  infinitesimal  fre- 
quency band  of  the  collision  broadened  line  may  be  viewed  as  being  broadened  by 

the  doppler  distribution.  The  resultant  lineshape,  referred  to  as  a Voigt  profile, 

45 

is  defined  by  Mitchell  and  Zemansky,  as 


exp  (-y^) 
a^  + (5  - y)^ 


dy 


(109) 


where  k is  defined  by  Eq.  (96),  5 
the  two  halfwidths. 


is  defined  by  Eq.  (97),  and  a is  the  ratio  of 


a 


tn  2 


and 


y = 


\/ln  2 , . 

— C (l/'  • V, 


(110) 


(111) 


Integration  over  the  dummy  frequency,  i/',  is  performed  in  order  to  include 
contributions  from  all  doppler-shifted  molecular  transitions  in  the  total  absorption 
coefficient  at  iv  - ).  For  the  present,  experimental  conditions  a = 0.22.3  - 0.  273. 

The  lineshapes  of  a [Jllire  doppler  transition  and  a Voigt  profile  with  a = 0.  25  are 
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plotted  in  Figure  43  normalized  to  the  peak  of  the  doppler  curve.  The  enhanced 
emission  in  the  wings  of  the  Voigt  profile  may  result  in  a significantly  different 
fractional  transmission  than  for  the  doppler  case.  The  fractional  transmission 
of  radiation  from  an  optically  thin  Voigt  broadened  source  attenuated  by  nonover- 
lapping Voigt  absorbers  is  given  by 


/ 
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L-00 


+ (f  - y) 


exp  I -k 


Np  L 4 
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+ (?  - y)^ 


dy 


(112) 


dr 


A solution  for  Eq.  (112)  to  permit  determination  of  t as  a function  L is  cur- 
rently being  developed.  The  function  r,  as  determined  from  Eq.  (112),  will  lie 
between  the  pure  doppler  and  collision  transmissions.  Equations  (102)  and  (108) 
have  been  used  to  evaluate  the  transmittance  of  the  C02(i'2)  band  in  the  doppler  and 
I.orentz  line  limits.  The  transition  linestrengths,  S,  and  frequencies  were  calcu- 
lated for  all  rotational  lines  of  the  fundamental  asymmetric  stretch  of  the 
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Figure  43.  Voigt  Llneshapes  for  a = 0 and 
a = 0.  23 
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O C O (626)  molecule  as  well  as  for  the  i/^  hot-band  transitions,  which 
involves  molecules  with  thermally  populated  excited  vibrational  modes,  and  for  all 
significant  bands  of  the  C ©2  , O C O , and  O C O isotopes.  The 
integrated  band  intensities,  S^,  of  the  transitions  included  in  the  calculation  are 
listed  in  Table  8 along  with  band  center  frequencies, 


Table  8.  €02(1/3)  Spectral  Properties 


1 

Isotope 

Transition 

Ground 

State 

^1^2  ^3 

s * 

V 

cm/molec. 

"^band* 

-1 

cm 

Branches 

Rotational  Levels 

J 

626 

00°0 

9.  598E-17 

2349. 146 

P,  R 

even  J only 

626 

01^0 

7.  66E-18 

2336.  637 

P,  Q.  R 

all  J,  c odd 

636 

00°0 

9.  602E-19 

2283, 490 

P,  R 

even  J 

628 

00°0 

3.33E-19 

2332.  112 

P,  R 

all  J 

626 

02^0 

3.08E-19 

2324. 148 

P,  Q,  R 

all  J,  c even 

626 

02°0 

1.  934E-19 

2327.  432 

P,R 

even  J 

626 

10°0 

1.  183E-19 

2326. 594 

P,  R 

even  J 

636 

01^0 

8.  18E-20 

2271.  763 

P,  Q,  R 

all  J,  c odd 

627 

00°0 

6.  37E-20 

2340. 01 

P,  R 

all  J 

628 

01^0 

2.  58E-20 

2324. 148 

P,  Q,  R 

all  J,  cSid  each 

626 

03^0 

1.  22E-20 

2311.  675 

P,  Q,  R 

all  J,  c odd 

627 

01^0 

4.  97E-21 

2327. 575 

P,  Q,  R 

all  J,  c&d  each 

636 

02^0 

3.46E-21 

2260.045 

P,  Q.  R 

all  J,  c even 



From  ref,  46 


All  transitions  observed  in  this  spectral  region  involve  A;  3 1,  and  may  have 

three  branches  labelled  P,  Q,  and  R for  rotational  quantum  changes  of  -1,  0,  and 
+ 1.  Molecules  with  a center  of  symmetry  (626,  6:i())  have  alternate  rotational  lines 
missing  from  all  transitions.  When  no  orbital  angular  momentum  exists  about  the 
internuclear  axis,  as  is  the  case  for  the  fundamental  transitions  of  all  the  isotopes, 
there  is  no  Q branch  evident  in  the  spectra.  When  angular  momentum,  f , does 


exist,  each  rotational  level  is  split  because  of  the  coupling  of  rotation  with  the 
orbital  angular  momentum.  The  nearly  degenerate  components  of  the  rotational 
level  are  labelled  c and  d. 

The  transition  line  center  frequencies  (in  cm“^)  involving  each  rotational  level 
(Including  each  of  the  nondegenerate  parity  doublets)  are  calculated  relative  to  the 
band  center  frequency  using  the  rotational  constants  given  in  McClatchey  et  al. 
The  total  integrated  bandstrength,  Sy,  is  then  divided  among  the  allowed  transi- 
tions, with  the  individual  lines trengths,  S,  given  by  McClatchey  et  al,  as 


Q 


rot 


exp  - 


S J'U'  + 1)  I 
( ^rot  ) ^band 


(113) 


is  the  rotational  partition  function  and  J'  is  the  initial  (upper)  rotational  level 

from  which  the  transition  originates.  S,  are  the  linestrength  factors  given  by  the 

37 

HOnl-London  Formulae  (Herzberg  ) as 
Sj^  = (J'  + 1)(J'  - f)/J' 

= (2J'  + 1)  f^/J'(J'  + 1)  <114) 

Sj^  = (J'  + 1 + f )(J'  + 1 - /)/(J'  + 1) 

which  are  very  good  approximations  to  the  true  rotational  linestrength  factors  for 
this  case  where  coupling  between  / and  J is  weak. 

These  linestrengths  are  then  used  to  calculate  the  fractional  transmission  of 
each  rotational  line  over  the  range  of  experimental  COg  pressures.  In  addition, 
the  emitted  intensity  per  molecule,  I,  for  each  rotational  line  is  calculated  by  inte- 
grating Eq.  (98) 


00  J 00 

-00  -00 

(115) 


' IbB  S(W/molecule-sr) 


46.  McClatchey,  R.A.,  Benedict,  W.  S. , Clough,  S.  A.,  Burch.  D.  E. , Calfee, 
R.F. , Fox,  K. , Rothman,  L.  S,  and  Caring,  J.S.  (1973)  AFCRL 
Atmospheric  Absorption  Line  Parameters  Compilation,  A.F.  S.  C.  Report 
AFCRL-TR-73-00^6. 
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The  total  emitted  intensity  at  each  of  the  experimental  CO2  pressures  was  obtained 
as  the  sum  of  the  values  for  each  rotational  line  for  both  the  pure  doppler  and 
collision  broadened  cases.  The  total  transmitted  intensity  per  molecule  for  all 
lines  is  plotted  in  Figure  44  as  a function  of  CO2  pressure  for  both  approxima- 
tions. (Note  the  predicted  intensity  is  radiation  per  CO2  molecule  at  a vibrational 
temperature  of  300*^K. ) Because  the  Lorentzian  .'ines  spread  more  of  the  total 
linestrength  into  the  wings  of  the  line,  transmission  by  the  collision  broadened 
lines  is  greater  than  for  the  doppler  case  for  all  CO2  pressures,  differing  by  a 
factor  of  five  at  1 torr.  In  the  limit  of  no  attenuation  (t  = D 


I 


^ ^ij  ’’ij  ^ ^%and.^ 

i i 3 •'  •' 

where  is  is  over  all  rotational  lines  and  j is  over  all  transitions.  This  equivalence 
is  maintained  to  within  2 percent  for  both  doppler  and  collision  broadened  cases. 


COp  Pressdre.  Torr 

Figure  44.  Total  Transmitted  Intensity  as  a 
Function  of  CO2  Pressure  for  Doppler  and 

Lorentz  Broadened  Lines.  P,  =9  torr 

Ar 
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The  role  of  the  hot  bands  and  isotopes  in  the  total  transmitted  intensity  is 
clarified  by  Figures  45  and  46.  In  Figure  45,  the  fractional  doppler  transmission, 

^ T of  each  of  the  six  strongest  transitions  in  Table  8 is  plotted  as  a function 
i 

of  CO2  pressure.  Doppler  radiation  emitted  by  the  626  fundamental  will  be 
severely  attenuated,  so  much  so  that,  even  though  its  integrated  band  strength  is 
over  twelve  times  greater  than  the  hot  band  (01^0)  Sy  value,  the  hot  band  radiation 
dominates  the  transmitted  doppler  radiation  spectrum  for  all  experimental  CO, 
pressures  as  shown  in  Figure  46.  At  pressures  near  0.  1 torr  of  COg,  the  01  0 
transition  of  the  626  isotope  will  account  for  nearly  half  to  the  total  transmitted 
intensify.  At  pressures  above  0.8  torr  even  the  636  fundamental,  whose  Sy  is  two 
orders  of  magnitude  less,  will  contribute  more  intensity  to  the  transmitted  spec- 
trum than  the  626  fundamental.  This  behavior  is  expected  when  the  exponential 
factor  dominates  Eq.  (99).  Also  plotted  in  Figure  45  is  the  fractional  transmission 
of  the  Lorentz  broadened  626  fundamental.  Self  absorption  for  the  Lorentz  case 
is  not  nearly  as  severe  and  the  626  fundamental  will  dominate  the  transmitted  rad- 
iation spectrum  over  the  COg  pressure  range  shown.  These  predictions  are  in 
good  agreement  with  published  values  of  the  626  fundamental  transmission  as  shown 
in  Table  9 (Bulos  et  al,  and  Kumer  et  aP°). 


Table  9.  Fractional  Transmission  of  626  Fundamental 


p 

C02 

0.  1 torr 

1.  0 torr 

1 

Doppler  i.imit  - this  work 

3.  24  X 10’^ 

47  1 

- Bulos  et  al 

3.  09  X 10’^ 

-- 

- Kumer  et  al^**  (225‘’K) 

3.  5 X 10'^ 

Lorentz  Limit  - this  work 

-- 

6.  5 X 10 

47 

- Bulos  et  al 

-- 

6.  1 X lO'^ 

47.  Bulos,  B.  H,,  and  F’helps,  A.  V.  (1976)  l^hys.  Rev.  A 14:615. 

48.  Kumer,  J.B. , and  James,  T.C.  (1974)  J.  Geophys.  lies.  79:638. 


Figure  45.  Fraction  of  Total  i nattenuated 
Band  Intensity  Transmitted  for  the  Six 
Strongest  Transitions  for  Doppler  Lines, 
and  for  the  Lorentz  Broadened  626  Funda- 
mental as  a Function  of  CO2  Pressure 
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Figure  46.  Relative  Contributions  of  Indi- 
vidual Bands  to  the  Total  Transmitted  Rad 
iation  as  a Function  of  C()2  Pressure 


CO.,  rrosstirf.  lorr 


It  is  to  be  noted  that  the  predictions  shown  in  Figure  46  may  be  used  to  justify 
the  neglect  of  absorption  by  ambient  CO.^  in  the  modeling.  This  absorption  will 
dominantly  effect  the  radiation  from  the  626  fundamental  band,  which  has  the 
strongest  bandstrength.  However,  as  shown  in  Figure  46  this  band  provides  only 
a small  fraction  of  the  total  transmitted  €02(1/^)  radiation  and  thus  a decrease  in 
its  transmissivity  will  not  significantly  effect  the  predicted  radiation 

levels.  Carrying  this  point  one  step  further,  although  it  is  expected  that  the  use 
of  a Voigt  profile  will  increase  the  transmission  of  the  626  fundamental  band, 
absorption  by  ambient  CO2  will  act  in  the  opposite  direction.  Thus  it  may  be  that 
the  use  of  a doppler  profile  provides  a realistic  estimate  of  the  total  transmitted 
radiation.  The  validity  of  this  conjecture  should  be  checked  by  more  detailed 
modeling. 

In  order  to  permit  comparison  with  the  experimental  spectra,  the  predicted 
transmitted  line  intensities  must  be  convolved  with  the  previously  described  sine 
slit  function  of  10  cm"'^  resolution  (see  Section  3.2).  During  convolution,  the  sine 
slit  function  was  truncated  at  the  second  zero  crossing  on  each  side  of  the  center 
frequency.  All  the  were  sorted  by  frequency,  summed  in  0.  4 cm  ^ cells,  and 
convolved  with  the  sine  function  to  determine  spectral  intensity  as  a function  of 
wavelength. 

Experimental  C02(i-'3)  spectra,  corrected  for  atmospheric  transmission  and 
detector  response,  are  displayed  in  Figures  47  and  48  for  nominal  CO2  pressures 
of  0.25  and  1 torn.  The  two  spectra  are  quite  similar,  exhibiting  a broad  double 
peak  in  the  region  of  maximum  intensity  and  two  secondary  peaks  at  = 2270  and 
2380  cm'^.  A small  bump  at  = 2360  cm'^  is  also  observed  in  the  data  for 
PCO2  1 torr.  Shown  for  comparison  in  Figures  49-53  are  the  predicted  spectra 
of  doppler  broadened  lines  for  P(202  ^ ^ respec- 

tively. The  prediction  of  Figure  49  corresponds  to  unattenuated  €02^’';^^  radiation 
and  can  be  seen  to  be  quite  different  from  the  data.  In  particular  the  spectral 
width  of  the  unattenuated  band  is  much  narrower,  the  double  peaks  in  the  region  of 
maximum  intensity  occur  at  different  wavenumbers  than  those  of  the  data,  the 
central  minimum  is  much  deeper,  and  no  secondary  peaks  arc  in  evidence.  It  is 
clear  that  the  effects  of  self  absorption  on  tlie  data  are  quite  severe.  In  contrast, 
both  data  and  prediction,  normalized  to  peak  intensity,  are  shown  in  Figures  51 
and  53  and  the  comparison  is  seen  to  be  excellent,  with  the  predictions  duplicating 
all  the  dominant  spectral  features  of  the  data. 

It  is  evident  from  comparison  of  Figure  49  to  Figures  50-53  that  the  effect  of 
increasing  CO2  pre.ssure  is  to  shift  the  spectral  intensity  to  the  red.  This  occurs 
because  the  hot  band  and  isotopic  band  radiation  increases  relative  to  that  of  the 
fundamental  band  (see  Figure  46).  The  predicted  contribution  of  each  of  the  domi- 
nant bands  (the  (>26  fundamental  band,  the  isotope  bands  and  the  hot  bands)  to  the 
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Figure  47.  Normalizeci  Experimental  Spectrum  at  t = 2.  6 msec  for 
Pco^  “ 25  torr,  I = 0.  9 mA,  V = 30  kV.  = 9 torr 
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Figure  48.  Normalized  Experimental  Spectrum  at  t 2.  fi  msec  for 
IV-/-,  ■ 1 torr,  I - 0.  3 mA,  v 30  kV,  Pa  „ - 9 torr 
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Figure  49.  Predicted  C02(u3)  Band  Intensity.  Doppler  broadened  line.s. 
Resolution  = 10  cm"^.  Pcc)2  ' ® torr,  no  attenuation.  Maximum 

intensity  - 4.39  X 10”^^  W/molecule  - cm  ^-sr 
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Figure  .'iO.  Predicted  Transmitted  C02(i'3)  Rand  Intensity.  Doppler  Broad 
ened  lines.  ^ ftesolution  10  cm"l.  Pc’().  tore.  Maximum  intensity 

2.  .591  y 10'2(>  VV' /molecule-sr-cm ^ 


"^200.00  2250-00  2300-00  2350-00  2400-00 

Wavenumber,  cm'l 

Figure  51.  Predicted  and  Measured  002(1/3)  Band  Transmitted  Intensity.  Doppler 
broadened  lines.  Resolution  = 10  cm'l.  PcOo  ^ torr.  Predicted  maximum 

intensity  = 1.  261  X 10*  W/molecule-sr-cm'l.  prediction, --data 
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Figure  53.  Predicted  and  Measured  C02(i'3)  Band  Transmitted  Intensity  of 
Doppler  Broadened  Lines.  All  transactions  included  with  self  absorption 
only.  Pc02  ’ torr.  Predicted  maximum  intensity  5.287  X 10~27 

W/molecule-sr-cm"^.  prediction,  data 


total  predicted  spectra  is  shown  in  Figures  54  and  55  for  COg  pressures  of  0.25 
and  1.  0 torr  respectively.  As  can  be  seen  the  626  fundamental  bands  provide  a 
dominant  contribution  only  to  the  secondary  peak  at  2380  cm  \ while  the  isotopic 
fundamental  bands  are  the  source  of  the  secondary  peak  at  2270  cm  Finally, 
the  primary  source  of  radiation  at  band  center  is  the  626  hot  bands.  The  fact  that 
all  these  spectral  features,  in  approximately  the  same  intensity  ratios,  are 
observed  in  the  data  to  be  taken  as  strong  evidence  that  the  present  modeling  pro- 
vides a good  approximation  to  the  actual  transmission  phenomena.  It  should  be 
noted  that  in  Figure  54  the  predicted  626  fundamental  intensity  at  = 2380  cm  ^ is 
slightly  larger  than  that  for  the  total  band  intensity.  This  is  an  artifact  introduced 
by  the  system  slit  function.  The  total  unconvolved  intensity  in  this  wavenumber 
region  is  of  course  slightly  greater  than  that  due  to  the  626  fundamental  alone; 
however,  the  effect  of  the  negative  side  lobes  in  the  system  slit  function  is  to  de- 
crease this  intensity  slightly  relative  to  that  of  the  626  fundamental. 
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Figure  54.  Comparison  of  Three  Sets  of  Doppler  Predictions  for  0021/3 
Band  Radiation  at  0.  25  torr  CO2  With  (a)  All  Transitions  of  Table  8 
Included  (Top  Trace),  (b)  No  Hot  Bands,  But  All  Isotopic  Fundamentals 
Included  (Middle  Trace),  and  (c)  Only  626  Isotope  Fundamental  Included 
(Lowest  Trace).  All  three  curves  are  normalized  to  peak  of  curve  a. 
Maximum  intensity  = 1.  261  X 10~26  w/molecule-sr-cm"l 


Finally  a set  of  spectral  transmission  predictions  for  Lorentz  broadened  lines, 
including  all  €02(1/2)  band  transitions  are  presented  in  Figures  56-59  for  CO2 
pressures  of  0.  1,  0.  25,  0.  5 and  1.  0 torr  respectively.  (The  predicted  unattenu- 
ated Lorentzian  €02(1/3)  spectrum  is  the  same  as  that  for  doppler  broadening  shown 
in  Figure  49).  It  can  be  seen  that  these  predicti^s  are  in  poor  agreement  with  the 
data  of  Figures  47  and  48,  thus  adding  further  credence  to  the  validity  of  the 
doppler  approximation. 

These  transmission  predictions  may  be  used  to  provide  an  estimate  of  the 
total  €02(1/3)  concentration  within  the  field  of  view.  The  total  observed  €02(13) 
radiation  1.^.  in  W/cm^-sr  may  be  related  to  the  €02(13)  concentration  by  the 
relationship 


N (hr) 
t'3 

^T  " 4irAfTRf 


(117) 
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Included  (Middle  Trace),  and  (c)  Only  62G  Isotope  Fundamental  Included 
(Lowest  Trace).  All  three  curves  are  normalized  to  peak  of  curve  a. 
Maximum  intensity  = 5.  287  X 10"27  W/molecule-sr-cm'l 


-20 

where  hr  is  the  vibrational  energy  of  the  rg  mode,  4.  6 X 10  J,  Tp  is  the  radia- 
tive lifetime  of  2,  3 msec,  Aj.  is  the  area  of  the  field  of  view,  52.  8 cm^,  and  f is  a 
correction  factor  for  transmission  losses  which  can  be  deduced  from  Figure  44. 

For  the  case  Pc02  ” torr,  I = 0.  9 mA,  which  was  analyzed  in  Section  3,  the 
actual  CO2  pressure  found  to  be  0.33  torr  and  thus,  from  Figure  44,  f 22.  As 
mentioned  earlier  the  data  is  not  absolutely  calibrated;  however,  an  order  of  mag- 
nitude estimate  of  the  COgtrg)  band  intensity  may  be  provided.  Using  this  estimate, 
Eq.  (117)  and  f = 22,  it  can  be  deduced  that 


N = 
‘'3 


1.  2 X lO^"*  molecules 


at  t 2.4  msec. 


(118) 


89 


Wavenumber,  cm'l 


Figure  56.  Predicted  Transmitted  Intensity  of  C02(i'3)  Radiation  for  Lorentz 
Broadened  Lines.  Pc02  ” ^Ar  = 9 torr.  Maximum  Intensity  = 

7.  689  X 10  W/molecule-sr-cm  ^ 
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Figure  57.  Predicted  Transmitted  Intensity  of  CC)2(i’3)  Radiation  for  Lorentz 
Broadened  Lines. 
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Figure  58.  Predicted  Transmitted  Intensity  of  CO^(i/3)  Radiation  for  Lorentz 
Broadened  Lines,  Rc02  " Torr.  Maximum  intensity  = 3.755  X 10"26 

W/molecule-sr-cm'l 
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Figure  59.  Predicted  Transmitted  Intensity  of  CC)2(i'3)  Hadiation  for  r.,orenf7. 
Broadened  I.ines.  Pcc)2  ' torr.  Maximum  intensity  2.703  X 10*2r) 

W / molecule -sr-c  III  ■ ^ 
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The  elobal  CO„(i'o)  excitation  rate  R.,  may  now  be  deduced  to  be  9.  2 X 10 

^23  i H .3 

mol/sec  from  Eqs.  (94),  (118)  and  the  observed  value  of  t ^ of  1.76  X 10  sec. 

Upon  comparing  this  with  the  total  ion  pair  creation  rate  within  the  field  of  view, 

as  determined  in  Eq.  (85),  it  can  be  seen  that  there  are  ~4.  9 002(1^2)  excitations 

per  ion  pair  created.  The  total  beam  powei  converted  to  002(1/2)  within  the  field 

of  view  is 

P -R  (hi/)  * 4.3  X 10"^  watts  , (119) 

"3  ''3 

as  compared  to  the  0.  8 W total  beam  power  deposited  within  the  field  of  view. 

Thus  ~ 5.  5 percent  of  the  beam  energy  deposited  in  the  gas  is  converted  to  002(1/2^ 
vibrational  energy.  Once  again  these  estimates  are  only  accurate  to  within  an 
order  of  magnitude. 

The  remaining  cases,  Pqq  ” ® ® been  studied  in 

detail.  However,  a crude  estimate  of  their  relative  peak  002(1/2^  concentrations 
may  be  developed  from  the  observed  peak  intensities.  Since  the  observed  C02(i'2) 
relaxation  times  are  approximately  the  same  between  0.  1-1.  torr,  the  C02(i'2) 
populations  will  be  roughly  proportional  to  the  peak  C02(v^)  intensity  multiplied 
by  the  predicted  correction  factor  f.  The  measured  peak  002(1/2^  intensity  and 
resulting  estimated  relative  002(1/2)  PoP'^^^tions  are  shown  in  Figure  60.  As  can 
be  seen,  the  estimated  002(1/2^  population  increases  only  some  50  percent  over  a 
factor  of  ten  change  in  CO2  pressure.  The  actual  variation  in  secondary  electron 
excitation  of  CO2  with  CO2  pressure  may  be  less  than  this,  inasmuch  as  the  effect 
of  vibrational  pumping  by  CO  has  not  been  taken  into  account. 

It  would  appear  then  that  the  percentage  of  deposited  beam  energy  which  is 
transferred  to  the  002(1/2)  state  does  not  vary  significantly  between  F’c'02  ' 

1.0  torr.  This  effect  is  most  probably  due  to  the  preponderance  of  argon  in  the 
gas  mix.  The  secondary  electrons  created  by  the  beam  will  typically  have  a sig- 
nificant amount  of  kinetic  energy  which  is  subsequently  lost  by  collisions  with  the 
gas  molecules.  Above  electron  energies  of  « 15  eV,  inelastic  electron-neutral 
collision  mechanisms  in  Ar  will  provide  the  dominant  electron  energy  decay.  How- 
ever, there  are  no  inelastic  channels  for  electron-argon  collisions  at  energies 
below  =!  15  eV  and  CO2  — even  strongly  diluted  in  A r — will  provide  the  dominant 
electron  energy  dissipation  for  this  energy  range.  Therefore,  it  is  possible  that 
a significantly  larger  percentage  of  beam  energy  might  be  channeled  to  the  C()2(i'2) 

state  at  higher  CO,  partial  pressures.  Indeed  it  has  been  estimated  by  Sawada 

49  ^ 

et  al  ' that  in  the  case  of  a thick  target  of  pure  CO2  ®20  percent  of  the  beam  energy 

49,  Sawada,  T.,  .Strickland,  D.J.,  and  Green,  A.E.S,  (1972)  J.  Geophvs.  Res. 

77:4812. 
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will  be  channeled  to  the  002(1'^)  mode.  It  should  be  stressed  that  the  estimates 
presented  in  Figure  60  are  crude  and  that  more  detailed  study  is  required. 


Figure  60.  Relative  Peak  €02(1^3)  Population  and 
Intensity  vs  CO2  Pressure.  P^r  ' ® torr. 

• - intensity,  O - population 


5.  SPVIVURV  \M)  CONCLUSIONS 

Both  002(1/3)  and  CO  fundamental  band  radiation  have  been  observed  in  elec- 
tron irradiated  mixtures  of  CO2/N2,  C02/Ar  and  C02/He.  In  all  cases,  the  evi- 
dent production  mechanism  for  the  vibrationally  excited  CO  is  the  electron/ion 
recombination  reaction 

C02^  + e - CO  + O . (26) 

The  CO  sp*  tra  for  one  case,  corresponding  to  beam  conditions  of  V = 30  kV, 

I = 0.9  mA,  and  gas  conditions  of  P^^,  = 9 torr,  Pco  ^ torr,  T - 300‘^K,  has 
been  reduced  in  some  detail.  The  rate  constants  for^he  reaction 
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CO(v)  + COg  - CO(v  - 1)  + 002(1/3) 


(3) 


have  been  determined  for  CO  levels  v = 1-G  and  lower  bounds  have  been  provided 

for  the  rate  constants  for  levels  v = 7-12.  The  rate  constants  for  reaction  (3), 

V 4-6,  as  determined  in  the  present  analysis,  are  a factor  of  2-3  larger  than 

24 

those  determined  by  Hancock  and  Smith.  The  cause  for  this  discrepancy  is  not 
clear. 

The  relative  creation  rates  for  CO(v)  production  during  beam  irradiation  have 
also  been  determined.  It  has  been  suggested  that  this  vibrationally  excited  CO  is 

most  probably  created  through  the  radiative  decay  or  quenching  of  excited  elec- 

34  42 

tronic  states  of  CO  which  are  known  ' to  be  produced  by  reaction  (26).  Although 
an  absolute  calibration  of  the  data  was  not  available,  it  has  been  crudely  (order  of 
magnitude)  estimated  that  one  vibrationally  excited  CO  molecule  is  created  per 
ion  pair.  Such  a production  efficiency  would  imply  that  ~3.  5 percent  of  the  elec- 
tron energy  transferred  to  the  gas  target  is  channeled  to  CO  vibrational  energy. 

The  C02(i-'3)  fluorescence  observed  in  the  Pc02  ’ torr  case  has  also  been 
examined.  The  dominant  002(1/3)  deactivation  mechanism  for  this  case  is 

002(13)  + Ar  - 003(31/2)  + Ar  (89) 

_ 1 c 3 

and  the  rate  constant  for  this  process  was  determined  to  be  1.  7 X 10  cm  /sec 

44 

in  agreement  with  previous  measurements  by  Yardley  and  Moore.  The  observed 
003(1/3)  radiation  was  determined  to  be  severly  optically  thick.  It  was  found  that 
the  transmitted  003(1/3)  spectra  could  be  well  approximated  by  an  absorption 
model  assuming  nonoverlapping  doppler  broadened  absorption  lines.  Using  this 
model  it  was  crudely  estimated  that  on  the  average  there  were  ~4.9  002(13) 
excitations  per  ion  pair  and  that  ~ 5.  5 percent  of  the  electron  energy  deposited  in 
the  gas  target  was  channeled  to  003(1/3)  vibrational  energy. 

It  must  be  stressed  that  these  results  are  specific  to  only  one  set  of  data, 
which  exhibited  rather  noisey  spectra,  and  should  be  confirmed  by  analysis  of 
additional  data  sets  taken  under  different  experimental  conditions.  A number  of 
approximate  techniques  used  in  the  analysis  should  also  be  examined  in  greater 
detail.  For  example,  the  effect  of  Voigt  line  shapes  on  the  002(1  3)  band  trans- 
mission function  should  be  determined  and  a technique  for  absolutely  calibrating 
the  data  should  be  developed. 
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